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Abstract 


New  methods  which  avoid  the  repeated  constructions  of  aggregate  Hamiltonian  matrices 
over  antisymmetric  basis  states  generally  required  in  conventional  calculations  of  adia¬ 
batic  potential  energy  surfaces  are  reported  for  ab  initio  studies  of  the  structures,  spectra, 
and  chemical  reactions  of  molecules  and  other  forms  of  matter.  A  representational  basis 
in  the  form  of  an  outer  spectral  product  of  atomic  eigenstates,  employed  in  the  absence 
of  overall  electron  antisymmetry,  is  shown  to  facilitate  development  of  an  exact  atomic- 
pair  expression  for  aggregate  Hamiltonian  matrices.  Unphysical  (non-Pauli)  eigenstates 
spanned  by  the  atomic  product  basis  are  identified  and  eliminated  by  a  unitary  trans¬ 
formation  of  the  Hamiltonian  matrix  obtained  from  the  matrix  representative  of  the  ag¬ 
gregate  electron  antisymmetrizer.  Hermitian  atomic  and  atomic-pair  interaction  matrices 
are  defined  which  individually  have  appropriate  asymptotic  separation  limits  and  can  be 
constructed  once  and  for  all  employing  unitary  transformations  of  antisymmetric  adiabatic 
diatomic  eigenstates  and  associated  potential  energy  curves.  The  aggregate  Hamiltonian 
matrix  constructed  in  this  way  includes  the  effects  of  overall  electron  antisymmetry  and 
incorporates  Wigner  rotation  matrices  for  representation  of  all  angular  dependencies.  A 
particular  implementation  of  the  theory  which  explicitly  enforces  the  limit  of  closure  in 
spectral-product  calculations  is  seen  to  correspond  to  adoption  of  canonically  orthogo- 
nalized  lineraly-independent  antisymmetrized  diatomic  states  obtained  from  conventional 
computational  procedures.  Illustrative  applications  to  the  well-known  low-lying  doublet 
and  quartet  states  of  the  H3  molecule  demonstrate  that  the  eigensurfaces  of  the  anti¬ 
symmetrizer  can  anticipate  the  structures  of  the  more  familiar  potential  energy  surfaces, 
which  are  found  to  be  in  satisfactory  accord  with  values  obtained  employing  prior  aggregate 
antisymmetrization.  Analyses  of  previously  devised  atoms-in-molecules  and  diatomic-in- 
molecules  methods  emphasize  the  importance  of  aggregate  basis-state  orthogonality,  use 
of  strictly  Hermitian  matrix  representatives,  and  unitary  transformations  of  atomic  and 
diatomic  fragment  information  in  such  connections,  attributes  which  have  generally  not 
been  included  in  applications  of  these  otherwise  promising  approaches.  The  present  de¬ 
velopment  is  seen  to  circumvent  the  limitations  of  these  early  semi-empirical  methods, 
providing  a  convergent  ab  initio  formalism  suitable  for  computational  applications. 
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1.0  Introduction 


Theoretical  studies  of  the  structures,  electromagnetic  spectra,  and  chemical  reactions  of 
molecules  and  other  forms  of  matter  require  accurate  approximations  to  the  complex  (Born- 
Oppenheimer)  potential  energy  surfaces  describing  their  ground  and  excited  electronic 
states  [1].  Such  surfaces  would  preferably  be  available  in  forms  that  allow  for  convenient 
solution  of  the  associated  molecular  Schrodinger  equation,  and  would  also  provide  for 
study  of  the  avoided  crossings  and  surface  intersections  that  can  give  rise  to  non-adiabatic 
dynamical  behaviors  [2].  The  ability  to  perform  simultaneous  quasi-classical  molecular 
dynamics  (MD)  and  electronic  structure  calculations,  in  which  the  electronic  surfaces  are 
determined  in  the  course  of  MD  simulation  at  the  particular  aggregate  configurations 
sampled  [3],  would  be  highly  desireable  in  cases  involving  sufficiently  heavy  atoms,  provided 
that  the  electronic  surfaces  employed  can  include  the  subtleties  associated  with  the  possible 
presence  of  large  numbers  of  interacting  aggregate  electron  configurations  [4] . 

Computational  methods  currently  employed  for  the  aforementioned  purposes  commonly 
entail  repeated  calculations  of  large  numbers  of  one-and  two-electron  integrals  over  explic¬ 
itly  antisymmetrized  basis  states  in  construction  of  a  many-electron  Hamiltonian  matrix 
for  a  range  of  atomic  spatial  arrangements,  followed  by  determination  of  the  energies  and 
associated  eigenfunctions  of  selected  electronic  states  at  each  atomic  arrangement  and 
construction  of  expectation  values  corresponding  to  physical  properties  of  interest.  Con¬ 
siderable  progress  has  been  made  in  implementing  and  applying  this  general  approach  in 
ab  initio  calculations  [5,6],  as  well  as  in  the  widely  employed  density-funtional  approxi¬ 
mation  [7].  Nevertheless,  it  would  be  highly  desireable  to  devise  new  methods  which  can 
avoid  the  repeated  calculations  required  in  the  conventional  approaches,  can  circumvent 
determinations  of  total  molecular  energies  and  their  associated  differencing  problems  in 
favor  of  calculations  of  constituent  atomic  energies  and  atomic-pair  interaction  energies, 
and  can  be  applied  to  both  ground  and  electronically-excited  states  on  a  common  basis. 

In  the  present  report,  an  atomic-pair-interaction  approach  to  ab-initio  studies  of  molecules 
and  other  aggregates  is  described  based  on  a  previously  described  atomic  spectral-product 
representation  of  the  electronic  degrees  of  freedom  [8-12],  The  development  accommodates 
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the  incorporation  of  ab-initio  diatomic  information  in  a  rigorous  fashion  free  of  the  subjec¬ 
tive  ad  hoc  assumptions  and  approximations  which  have  rendered  such  fragment-based  ap¬ 
proaches  largely  semi-empirical  in  nature  [13,14].  The  present  method  is  ultimately  based 
on  conventional  variational  calculations  in  L 2  representations  of  many-electron  states,  guar¬ 
anteeing  its  bounded  (Hyllcraas-Undheim)  convergence  when  closure  is  achieved,  in  which 
limit  the  distinction  between  use  of  atomic-product  and  explicitly  antisymmetrized  aggre¬ 
gate  basis  states  is  shown  to  become  inconsequential.  As  presented  here,  the  theory  is 
applicable  to  finite  atomic  aggregates  which  dissociate  into  neutral  species  on  their  lowest- 
lying  potential  energy  surfaces,  but  which  can  involve  arbitrary  admixtures  of  covalent, 
ionic,  van  der  Waals,  and  metallic  interatomic  interactions.  Applications  of  the  approach 
to  ionic  systems,  or  to  neutral  systems  which  can  give  rise  to  ion-pair  states  in  asymp¬ 
totic  separation  limits,  require  only  minor,  largely  self-evident,  modifications  of  the  present 
development. 

The  atomic  spectral-product  basis  employed  in  the  development  is  known  to  span  the  to¬ 
tally  antisymmetric  irreducible  representation  of  the  aggregate  electron  symmetric  group 
once  and  only  once,  making  it  suitable  for  analytical  and  computational  studies  of  the 
electronic  structures  of  molecules  and  other  forms  of  matter  [8-12],  The  matrix  represen¬ 
tative  of  the  aggregate  electron  antisymmetrizer  constructed  in  the  spectral-product  basis 
can  be  employed  in  separating  the  totally  antisymmetric  and  non-totally-antisymmetric 
(non-Pauli)  states  spanned  by  the  basis,  and  in  correspondingly  isolating  the  physically 
significant  block  of  the  Hamiltonian  matrix  by  unitary  transformation.  Hermitian  matrix 
representations  of  atomic  and  atomic-pair  operators  are  devised  in  this  way  which  individ¬ 
ually  have  well-defined  asymptotic  separation  limits,  and  which  need  be  constructed  only 
once  for  multiple  applications.  These  attributes  of  the  approach  facilitate  development  of 
an  ab  initio  unitary  transformation  formalism  which  provides  an  exact  atomic-pair  rep¬ 
resentation  of  aggregate  Hamiltonian  matrices,  including  particularly  modifications  due 
to  aggregate  electron  antisymmetry  of  the  bare  atomic-pair  interactions  consequent  of 
their  incorporation  in  the  atomic  aggregate.  The  development  furthermore  allows  for  the 
incorporation  of  Wigner  rotation  matrices  in  an  efficient  representation  of  all  angular  de¬ 
pendences  of  the  interaction  matrix  elements  [15]. 
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A  particular  computational  implementation  of  the  theoretical  development  is  described 
which  enforces  closure  limits  over  the  formally  infinite  spectral  sums  required  to  achieve 
wave  function  antisymmetry,  providing  expressions  which  can  be  evaluated  with  largely 
standard  ab-initio  computational  methods  [5-7].  Illustrative  applications  in  the  case  of  the 
low-lying  doublet  and  quartet  states  of  the  H3  molecule  demonstrate  that  the  eigensurfaces 
of  the  electron  antisymmetrizer  can  anticipate  the  structures  of  the  corresponding  more 
familiar  potential  energy  surfaces,  which  are  found  to  be  in  accord  with  values  obtained 
employing  conventional  computational  methods.  This  particular  implementation  of  the 
theory  also  helps  to  identify  the  origins  of  complications  encountered  in  applications  of 
related  semi-empirical  atoms-in-molecules  [13]  and  diatomic-in-molecules  [14]  formalism 
based  on  explicitly  antisymmetric  aggregate  basis  states,  which  issues  have  apparently 
precluded  their  widespread  adoption  to  date  as  viable  computational  approaches.  Modifi¬ 
cations  of  these  approaches  based  on  the  present  analysis  are  seen  to  lead  to  potentially 
convergent  ab  initio  computational  formalisms 

The  theory  is  described  in  Section  2,  where  the  atomic  spectral  representation  of  aggregate 
electronic  degrees  of  freedom  is  defined,  the  unitary  transformation  formalism  for  isola¬ 
tion  of  the  physical  portion  of  the  spectrum  and  for  incorporation  of  ab  initio  diatomic 
information  is  described,  and  the  exact  atomic-pair  form  of  the  aggregate  Hamiltonian 
matrix  reported.  The  particular  computational  implementation  of  the  approach  devised 
is  described  and  illustrative  calculations  of  the  doublet  and  quartet  potential  energy  sur¬ 
faces  of  H3  provided  in  Section  3.  Comparative  and  concluding  remarks  made  in  Section  4 
place  the  present  development  in  the  context  of  the  aforementioned  earlier  semi-empirical 
approaches  [13,14],  help  to  clarify  the  origins  of  complications  encountered  in  their  appli¬ 
cations,  and  demonstrate  how  the  present  development  circumvents  these  difficulties. 

2.0  Theory 

The  atomic  spectral-product  approach  to  aggregate  electronic  structure  is  described  in 
Section  2.1,  the  unitary  transformation  formalism  for  incorporation  of  diatomic  information 
is  presented  in  Section  2.2,  and  the  resulting  exact  atomic-pair  form  of  the  aggregate 
Hamiltonian  matrix  is  reported  in  Section  2.3. 
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2.1  Atomic  Spectral-Product  Formalism 

Solutions  of  the  Schrodinger  equation  [16] 

H(r  :  R)^(r  :  R)  =  ^(r  :  R)  -E{R)  (1) 


for  an  IV-atom  molecule  or  other  finite  atomic  aggregate  with  Hamiltonian  operator 

N  N 

H(r  :  R)  =  J2  {«■'">(*)  +  J2  V{a’0){i-,j-Rct 3)}  (2) 

a=l  /3=a+ 1 

can  be  constructed  employing  a  formally  complete  orthonormal  square-integrable  ( L 2 ) 
outer-spectral-product  representation  of  the  form  [8-12] 


<F(r  :  R)  =  {$(1)(i)  <g>  $(2)(2)  <g>  ■  •  •  $(A°(n)} 


The  Hamiltonian  operator  H(r  :  R)  in  Eqs.  (1)  and  (2)  includes  a  sum  over  atomic  terms 


2 


e2 


+  £  h 


i'=i+ 1 


and  a  sum  over  pairwise- atomic  Coulombic  interaction  terms  [17] 


7  7  „2  na  7  2  n/3  7  2  n«  n/3  2 

:  **)  =  h&L.  -  £  ^  -  X  ^  +  ££  h 


where  r  and  R  refer  to  all  nt  aggregate  electron  and  N  atomic-position  coordinates,  re¬ 
spectively.  The  desired  row  vector  of  eigenfunctions  T'f  r  :  R)  and  the  diagonal  matrix  of 
associated  eigenvalues  E (R)  of  Eq.  (1)  generally  include  both  physically  significant  and 
unphysical  (non-Pauli)  solutions  [12].  The  row  vectors  in  Eq.  (3)  comprise  the 

orthonormal  antisymmetric  multiplet  eigenstates  of  the  constituent  atoms  (a:  =  i  to  N ) 
specified  by  the  quantum  numbers  (E,  L,  Ml,  S,  Ms,  P)a  [18],  with  i  (j)  referring  to  the 
coordinates  relative  to  the  atomic  center  Ra  {Rp)  of  the  na  (np)  electrons  arbitrarily 
assigned  to  the  atom  a  (/?).  The  atomic  eigenstates  can  be  variationally  determined  in 
L2  basis-state  representations  and  need  not  be  exact  eigenfunctions  of  the  corresponding 
atomic  Hamiltonian  operators  for  the  ensuing  development  to  converge,  so  long  as  the  de- 
numerably  infinite  outer-product  (®)  Hilbert  space  of  Eq.  (3)  includes  a  suitable  closure 
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limit  over  the  range  of  atomic  configurations  R  of  interest  [19].  The  subscript  “O”  in 
Eq.  (3)  indicates  the  adoption  of  a  particular  ordering  convention  for  the  vector  sequence 
of  spectral-product  functions  employed,  the  consequences  of  which  convention  are  indi¬ 
cated  when  appropriate  in  the  sequel.  More  detailed  descriptions  of  this  and  of  the  other 
notational  conventions  of  Eqs.  (1)  to  (5)  are  reported  elsewhere  [8-12], 

Employing  the  formally  complete  orthonormal  basis  of  Eq.  (3)  in  variational  solution  of 
Eq.  (1)  gives  the  familiar  matrix  Schrodinger  equation 

H (R)  •  UH (R)  =  UH(Ji)  •  E (R),  (6) 

where  the  columns  of  the  unitary  matrix  Uh(-R)  which  diagonalizes  the  Hermitian  Hamil¬ 
tonian  matrix  H(il)  provide  the  desired  eigenvectors  \I/(r  :  R)  =  $(r  :  R )  •  Uh(-R) 
following  commonly  employed  labelling  and  ordering  conventions  [16].  The  Hamiltonian 
matrix  of  Eq.  (6), 

N  N 

■H(R)  =  (<t>(r:R)\H(r:R)\<t>(r:R))  =  J2{H{a)+  E  (7) 

a=l  /3=a+ 1 

includes  constant  atomic  terms  =  (<I>(r  :  R)\H^a',(i)\^(r  :  R))  and  atomic-pair 

interaction  terms  \^a^\Rap)  =  ($(r  :  R)\V('a,l3\r1  j  :  Rap)\&(r  :  R))  which  are  the 
Hermitian  matrix  representatives  of  the  corresponding  atomic  [Eq.  (4)]  and  interaction  [Eq. 
(5)]  operators  in  the  aggregate  Hamiltonian  operator  of  Eq.  (2).  The  atomic  energy  matrix 
includes  only  diagonal  multiplet  energies  which  are  independent  of  atomic  position, 
consequent  of  choosing  the  orthonormal  atomic  states  <I>i:rd  (i)  as  eigenstates  or  pseudo¬ 
eigenstates  of  the  respective  atomic  Hamiltonian  operators,  whereas  the  interaction  matrix 
is  generally  non-diagonal  and  depends  explicitly  upon  the  vector  separation 
Rap  of  the  two  indicated  atoms,  but  not  upon  the  individual  laboratory-frame  positions  of 
the  two  atoms,  nor  upon  the  position  vectors  of  the  other  (N  —  2)  atoms  in  the  aggregate. 

The  particularly  simple  form  of  Eq.  (7)  is  largely  a  consequence  of  the  orthonormality  of 
the  spectral-product  basis  employed  in  the  absence  of  prior  enforcement  of  overall  aggregate 
electron  antisymmetry,  the  use  of  atomic  multiplet  eigenstates  in  the  representation,  and 
the  atomic  pairwise-additive  nature  of  the  interaction  terms  in  the  Hamiltonian  operator 
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of  Eq.  (2).  It  lias  been  shown  that  the  basis  of  Eq.  (3)  spans  the  totally  antisymmetric 
representation  of  the  aggregate  n /-electron  permutation  (or  symmetric  -  Snt)  group  once 
and  only  once  in  the  limit  of  closure,  although  other  unphysical  non-totally-antisymmetric 
(non-Pauli)  representations  of  Snt  are  also  spanned  by  the  basis  [12,20-24], 

The  totally  antisymmetric  eigenstates  and  the  non-Pauli  eigenstates  obtained  from  Eqs. 
(1)  to  (7)  can  be  separated  by  constructing  the  transformed  Hamiltonian  matrix 

N  N 

Hs(fl)  =  UsM*  ■  H(fl)  •  Us(fl)  =  J2  {  H^(«)  +  J2  Vs 

a=l  /3=a+l 

__  f  {Hs(-R)}pP  {Rs(R)}pu\  f{Us(R)}pp  0  \ 

V(HS  (R)}up  {Us(R)}nuJ  V  0  m 

where  H ^ (R)  =  Us(JR)t  •  H(Q)  •  US(R)  and  V^’0)(R)  =  US{R)]  •  V^\Raf3)  ■  US(R) 
refer  to  transformations  of  the  individual  atomic  and  interaction  matrices  in  Eq.  (7), 
{Hs(H)}pp  and  {Hs(-R)}wu  are  the  physical  (pp)  and  unphysical  (uu)  blocks  of  the  trans¬ 
formed  Hamiltonian  matrix  which  provide  the  antisymmetric  and  the  non-Pauli  solutions, 
respectively,  and  the  off-diagonal  blocks  vanish  in  the  limit  (— »)  of  spectral  closure  [9-12], 

The  unitary  transformation  matrix  Ug(-R)  required  to  accomplish  the  partitioning  of  Eq. 
(8)  is  obtained  from  diagonalization 

Us,fl)t.S(fl).Us(K)  =  Sd(K)=({S^»»  {Sd{%uy(Q10"  o)  <9> 

of  the  Hermitian  matrix  representative  of  the  antisymmetrizer  P/i  constructed  in  the 
spectral-product  basis, 

S (R)  =  Q1/2($(r  :  ii)|PA|^(r  :  R))  =  (PA^(r  :  R)\PA$(r  :  R)),  (10) 

which  is  seen  to  be  proportional  to  the  metric  matrix  of  the  explicitly  antisymmetrized 
spectral-product  basis.  The  normalization  of  the  antisymmetrizer  employed  in  the  present 
development,  PA  =  (rp!)-1/2(ni!n2!  •  •  -nAr!)-1^2  Sp=i(— P5pPpi  is  chosen  to  insure  that 
S  {R  — >  oo)  — >  I  when  all  the  atomic  states  of  Eq.  (3)  from  the  different  atoms  no  longer 
“overlap”  in  the  sense  of  Eq.  (10)  in  the  complete  dissociation  limit.  Adoption  of  this 
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normalization  convention,  which  is  commonly  employed  to  insure  unity  normalization  of 
the  antisymmetrized  product  basis  :  R)  in  the  limit  ( R  — >  oo)  [13],  gives  the  value 

Q  =  nt\/(ni\ri2\  ■  ■  -njv!)  for  the  redundancy  factor  of  the  antisymmetrized  basis  appearing 
in  Eq.  (10)  [22], 

The  non- negative  eigenvalue  matrix  S d(R)  in  Eq-  (9)  is  partitioned  into  an  upper  diago¬ 
nal  {Sd{R)}pP  containing  the  largest  eigenvalues  ordered  in  decreasing  value,  which  tend 
to  the  indicated  common  redundancy  factor  Q  in  the  closure  limit,  and  a  lower  diagonal 
{S d{R)}uu  °f  eigenvalues  which  tend  to  zero  in  this  limit  [9-12],  The  eigenvalues  are  so 
ordered  by  appropriate  arrangement  of  the  columns  of  the  transformation  matrix  Us(-R), 
whereas  the  dimensions  of  the  blocks  {Hs(-R)}pp/MU  of  Eq.  (8)  in  finite  order  are  deter¬ 
mined  by  that  block  of  the  diagonal  matrix  Sa(R)  of  Eq.  (9)  subjectively  selected  as  its 
physical  part  in  the  particular  atomic  spectral-product  representation  employed. 

The  unitary  transformation  matrix  Us(-R)  of  Eq.  (9)  obtained  from  the  metric  matrix 
of  Eq.  (10)  is  responsible  for  incorporating  the  non-local  effects  of  inter-atomic  aggregate 
electron  permutation  symmetry  in  the  Hamiltonian  matrix  of  Eq.  (7)  by  virture  of  Eq.  (8) 
[9-12],  That  is,  the  diagonalization  of  Eq.  (9)  corresponds  to  construction  of  eigenstates 
of  the  Hermitian  antisymmetrizer  in  the  spectral-product  basis,  so  that  the  columns  of  the 
transformation  matrix  Us(-R)  provide  the  eigenstates  of  the  operator  Pa  in  the  spectral- 
product  representation  in  the  form 

&s(r  :  R)  =  Hr  :  R)  ■  U S{R)  -  {{$s(r  :  R)}p,  {H (r  :  R)}u}-  (11) 

Here,  { <I>g(r  :  R)}p  contains  the  totally  antisymmetric  states  corresponding  to  the  non-zero 
eigenvalues  of  the  antisymmetrizer  and  { s ( r  ■  R)}u  contains  the  non-Pauli  states  corre¬ 
sponding  to  the  zero  eigenvalues  of  the  antisymmetrizer  in  the  closure  limit.  Consequently, 
it  follows  that  Eq.  (8)  constitutes  a  unitary  transformation  of  the  Hamiltonian  matrix 
of  Eq.  (7)  constructed  in  the  spectral-product  basis  to  the  orthonormal  permutation- 
symmetry-adapted  representation  of  Eq.  (11). 

Equations  (9)  to  (11)  indicate  that  the  spectrum  of  Pa  acting  in  the  domain  of  the  spectral- 
product  basis  of  Eq.  (3)  is  that  of  a  compact  Hermitian  operator  [25],  with  zero  providing 
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a  lower  limiting  point  of  accumulation  of  the  spectrum  and  the  associated  non-Pauli  states 
{‘I’sty  :  i?)}M  of  Eq.  (11)  corrersponding  to  its  null  space  Pn{^ s{r  '■  R)}u  0,  whereas 
the  antisymmetric  states  {3>s(^  :  R)}p  of  Eq.  (11)  are  associated  with  the  upper  limiting 
point,  Pa{®s (r  :  R)}p  -►  Q1/2{$s(r-  :  R)}p,  and  correspond  to  its  non-vanishing  point 
spectrum. 

The  condition  S (R  — >  cx))  — »  I  indicated  above  can  always  be  satisfied  for  a  sufficiently 
large  dissociation  limit  (R  — >  oo)  which  separates  all  atoms  beyond  the  range  of  all  atomic 
eigenstate  overlap  in  any  particular  denumerable  L 2  spectral-product  representation.  In 
this  limit,  the  transformation  matrix  of  Eq.  (9)  becomes  a  constant  unitary  matrix  Us(oo) 
[=  Us  (R  —>  oo)],  the  physical  portion  of  Eq.  (11)  becomes  a  unitary  transformation  of 
the  entire  spectral-product  basis  of  Eq.  (3)  [<J>s(r  :  R  — >  oo)  — >  <J?(r  :  oo)  •  Us(oo)],  the 
unphysical  portion  of  Eq.  (11)  is  correspondingly  null,  the  interaction  matrices  of  Eq.  (7) 
vanish  \y(a,l3\Rap  — >  oo)  -0],  and  the  Hamiltonian  matrix  of  Eq.  (8)  becomes 

N 

US(R  ->  oo)  ->  Us(oo)f  •  H  (R  ->  oo)  •  Us(oo)  =  Us(oo)f  •  {  H^j  •  Us(oo).  (12) 

a=l 

In  view  of  the  unitary  nature  of  the  transformation  of  Eq.  (12),  the  aggregate  energies 
in  the  dissociation  limit  are  seen  to  be  sums  of  the  individual  non-interacting  atomic 
energies,  as  anticipated.  Moreover,  although  the  unitary  matrix  Us(oo)  obtained  from 
Eq.  (9)  is  formally  arbitrary  in  this  limit,  as  a  consequence  of  the  degeneracy  of  the 
eigenvalues  of  the  metric  matrix  [S (R  -»■  oo)  ->  I],  it  is  possible  to  choose  the  aggregate 
eigenfunctions  <&(r  :  oo)  •  Us(oo)  as  states  of  good  total  aggregate  spin  multiplicity,  and 
of  other  possible  aggregate  point-group  symmetry,  insuring  that  Us(oo)  commutes  with 
the  indicated  atomic  matrix  sum  in  Eq.  (12).  In  this  case  the  Hamiltonian  matrix  of  Eq. 
(12)  is,  of  course,  identical  with  the  indicated  sum  of  unperturbed  atomic  energies. 

The  absence  of  explicit  inter-atomic  antisymmetry  in  the  spectral-product  representation 
of  Eq.  (3)  and  the  condition  S(R  — >  oo)  — »  I  insures  that  well-defined  Hermitian  repre¬ 
sentatives  Hg*-*  (R)  and  Vg “’^(i?)  of  the  atomic  and  interaction-energy  operators  of  Eqs. 
(4)  and  (5),  respectively,  are  obtained  which  individually  have  appropriate  asymptotic- 
separation  limits,  providing  the  correct  atomic-product  aggregate  eigenstates  and  the  as- 
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sociated  sums  of  atomic  energies  in  Eq.  (12)  in  this  limit.  By  contrast,  an  explicitly 
antisymmetrized  basis  P,\  d*  ( v  '■  R)  does  not  allow  for  the  definition  of  individual  atomic 
and  atomic-interaction  operators  and  associated  Hermitian  matrix  representatives  from 
the  expectation  value  of  the  aggregate  Hamiltonian  operator  which  individually  have  ap¬ 
propriate  aggregate  dissociation  limits,  the  various  potential  energy  terms  contributing 
differently  as  intra-  or  inter-atomic  interactions  consequent  of  the  different  atomic  assign¬ 
ments  of  electrons  arising  from  the  different  permutations  in  the  basis  states. 

Of  course,  the  total  Hamiltonian  matrix  in  the  antisymmetrized  basis  has  a  satisfactory  dis¬ 
sociation  limit  ( R  — >  oo),  although  this  matrix  is  generally  singular  in  the  spectral  closure 
limit  consequent  of  the  linear  dependence  of  the  overcomplete  explicitly  antisymmetrized 
basis  at  any  finite  R  [12,20,22],  Additionally,  the  total  or  aggregate  electron  probability 
density  matrices  <I>(r  :  Rp  ■  $(r  :  R)  and  PA*&{r  :  Rp  ■  pA*&{r  :  R)  are  identical  in 
the  limit  R  — >  oo,  insuring  the  distinguishability  of  the  electrons  assigned  to  the  different 
atoms  in  this  limit  on  basis  of  11011-overlap  of  the  atomic  states  employed,  even  in  the 
presence  of  explicit  antisymmetry  [26].  This  equivalence  by  itself,  however,  is  insufficient 
to  insure  that  the  explicitly  antisymmetric  representation  provides  atomic  and  interaction 
operators  and  matrices  which  individually  have  appropriate  asymptotic  separation  limits. 

2.2  Atomic-Pair  Transformation  Formalism 

The  particularly  simple  forms  of  the  atomic  and  interaction-energy  terms  appearing  in  the 
Hamiltonian  matrix  of  Eq.  (7)  suggest  the  possibility  of  systematic  evaluation  of  these 
component  matrices  in  a  series  of  individual  atomic  and  atomic-pair  calculations  which 
can  be  performed  once  and  for  all  independent  of  any  particular  atomic  aggregate  arrange¬ 
ment.  Since  atomic  and  diatomic  eigenspectra  are  generally  available  from  experimental 
and  previous  theoretical  studies,  the  quality  of  the  required  calculations  of  component  in¬ 
formation  can  be  critically  evaluated  in  this  approach  through  suitable  comparisons  prior 
to  their  incorporation  in  the  spectral-product  representation  of  Eqs.  (1)  to  (12). 

In  order  to  explore  the  foregoing  line  of  reasoning,  it  is  helpful  to  write  the  Hamiltonian 
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operator  of  Eq.  (2)  in  the  three  alternative  but  equivalent  forms  [14] 

N  N  (AT 

H(r  :R)  =  Y.  E  {£<“jS)(i. 3  ■  Raf>)  ~  ■  Rap  -  °o)} 

a=l  /3=a+ 1  ' 

N  N  /  qx 

=  E  E 

a=l  /3=cc+l  ' 

N  N 

=  E{^(“’(i)+  E  (13) 

a=l  /3=a+l 

where 

H{a’P\iJ  :  iia/j)  =  i7(a)(*)  +  HW(j)  +  :  iJa/3)  (14) 

is  the  diatomic  Hamiltonian  operator  for  atoms  a  and  (3.  The  operator  of  Eq.  (14)  is 
a  symmetric  function  of  electron  coordinates  i  and  j,  and,  accordingly,  has  a  physically 
significant  Hermitian  matrix  representative  in  the  spectral-product  basis,  whereas  the  in¬ 
teraction  operator  V(a,P\i;j  :  Ra/3 )  of  Eq.  (5)  is  not  a  symmetric  function  of  electron 
coordinates  and,  consequently,  has  no  physical  significance  [27]. 

Employing  Eq.  (13),  the  Hamiltonian  matrix  of  Eq.  (7)  can  be  written  in  the  three 
corresponding  equivalent  forms 

N  N  ,  , 

H(«)  =  E  E  {H<“,'5)(«^)-7773TTH<“-'5)(-R=/S^~)} 

a=l  /3=a+l 

N  N  ,  , 

=  E  E  HWIw-k^{HW+HW} 

a=l  /3=a+l  ' 

N  N 

=  £  {  H(Q°  +  ^  (15) 

a=l  /3=a+l 

where  the  terms  in  the  double  sum  in  the  last  line  define  the  (ct,  /3)  interaction  potential 
matrices  through  the  expression  [cf.,  Eq.  (14)] 

H ^\Rap)  =  ($(r  :  R)\H^\iJ  :  Rap)\$(r  :  R)) 

=  H(Q)  +H(/3)  +  V(a^\Ra/3).  (16) 
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Although  the  diatomic  matrix  of  Eq.  (16)  as  written  has  the  dimension  of  the  entire 
spectral-product  basis  of  Eq.  (3),  it  is,  of  course,  determined  by  the  diatomic  submatrix 
constructed  from  integrals  over  the  subspace  of  atomic-pair  functions  <!>(“)  (i)  (8)  in 

which  the  diatomic  operator  of  Eq.  (14)  acts,  and  by  the  outer  product  of  this  submatrix 
with  the  unit  matrices  expressing  the  orthonormality  of  the  remaining  IV— 2  atomic  spectral 
basis  states  of  Eq.  (3).  Moreover,  in  view  of  the  ordering  convention  of  Eq.  (3),  the 
matrices  tl^a’^\Rap)  for  the  different  atomic  pairs  are  each  differently  arranged  in  terms 
of  their  atomic  and  pair-interaction  matrix  elements  [8,12], 

The  transformed  Hamiltonian  matrix  of  Eq.  (8)  in  the  permutation-symmetry-adapted 
basis  of  Eq.  (11)  takes  the  three  equivalent  forms 

Hs(^)  =  E  £  {Hga,/3)(.R)  —  YyHga,/3)(i2  :  Rap  — >  oo)} 

a=l  /3=a+ 1 

N  N  rATr,\ 

TE  {4“'«(B)-  E-d{HW(R)  +  HW(H)}} 

a=l  /3=a+l 

N  N 

=  E{Hs“,(-R>+  E  {«?•«(«)  -H^(fl)-Hsm(S))},  (17) 

a= 1  f3=a-\-l 

where,  as  in  Eq.  (15),  the  terms  in  brackets  in  the  last  line  define  the  (a,  (3)  interaction 
matrices  \^’^\r)  through  the  expression  [cf.,  Eqs.  (8)  and  (16)] 

H ^\R)  =  Us(JR)t  ■  U^\Raf3)  •  Us (R) 

=  ($s (r  :  R)\H^\i,j  :  ^a/3)|$s(r  :  R)) 

-  H  {s\R)  +  H{s\R)  +  V(s,I3){R)  (18) 

for  the  corresponding  transformed  diatomic  Hamiltonian  matrix.  The  notation  employed 
in  the  second  term  in  the  first  lines  of  Eq.  (17)  means  that  the  limit  Rap  — >  oo  is  to 
be  taken  in  the  operator  :  Rap)  and  in  the  matrix  H (a,P\Rap),  but  not  in 

the  functions  $s(r>  :  R)  or  the  transformation  matrix  Us(-R)  in  Eq.  (18).  Accordingly, 
with  this  convention  the  second  term  is  seen  to  reduce  correctly  to  the  indicated  sum  of 
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transformed  atomic  Hamiltonian  matrices  required  to  cancel  the  additional  atomic  terms 
included  in  the  first  double  sum  over  atomic  pairs,  as  in  Eq.  (15). 

The  diatomic  Hamiltonian  matrix  of  Eq.  (18)  is  seen  to  include  modifications  of  the  corre¬ 
sponding  bare  diatomic  matrix  of  Eq.  (16)  due  to  the  inter-atomic  antisymmetry  required 
upon  immersion  of  the  interacting  atomic  pair  in  the  molecular  aggregate,  accomplished 
through  the  action  of  the  transformation  matrix  Us(-R)  of  Eq.  (9).  Accordingly,  Eqs.  (17) 
and  (18)  can  provide  an  exact  representation  of  the  entire  aggregate  Hamiltonian  matrix 
in  the  limit  of  closure  in  the  form  of  sums  of  Hermitian  atomic  and  atomic-pair-interaction 
matrices  which  individually  have  appropriate  dissociation  limits.  In  employing  these  ex¬ 
pressions  in  computational  applications,  it  is  convenient  to  first  perform  an  eigenvalue 
decomposition  of  the  bare  Hamiltonian  matrix  of  Eq.  (16)  and  to  retain  the  eigenstates 
and  energies  so  obtained  for  subsequent  aggregate  applications,  rather  than  to  retain  the 
diatomic  Hamiltonian  matrix  itself.  In  this  way,  as  will  be  seen  below,  conventional  compu¬ 
tational  methods  can  be  brought  to  bear  most  conveniently  in  evaluations  of  the  aggregate 
Hamiltonian  matrix  of  Eqs.  (15)  to  (18). 

Incorporation  of  diatomic  fragment  information  in  the  foregoing  exact  atomic-pair  ex¬ 
pressions  as  an  alternative  to  direct  evaluation  of  the  interaction  terms  in  the  aggregate 
Hamiltonian  matrix  of  Eqs.  (7)  and  (8)  is  complicated  by  a  number  of  issues.  First,  di¬ 
atomic  calculations  are  most  conveniently  performed  in  collinear  coordinate  systems  having 
their  z-axes,  for  example,  aligned  and  employed  as  a  common  axis  of  quantization,  rather 
than  in  the  commonly-oriented  coordinate  systems  describing  the  individual  atoms  a  and 
j3  in  the  aggregate  laboratory  coordinate  frame.  Second,  the  spectral-product  basis  of 
Eq.  (3)  will  generally  span  unphysical  states  of  the  diatomic  Hamiltonian  operator  of  Eq. 
(14),  suggesting  that  the  spectrum  of  the  Hamiltonian  matrix  of  Eq.  (16)  will  include 
both  physical  and  non-Pauli  eigensolutions,  both  of  which  can  potentially  contribute  to 
the  development.  Third,  the  transformations  expressing  the  aggregate  Hamiltonian  ma¬ 
trix  in  terms  of  its  component  diatomic  parts  require  construction  in  the  absence  of  any 
possible  ambiguity  or  inconsistency  which  might  preclude  applications  of  the  method  as  a 
convergent  ab  initio  computational  approach. 
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2.2.1  Spatial  Frame  Transformation 


To  accomplish  the  transformation  from  aggregate  to  diatomic  frame,  a  row  vector  of 
spectral-product  functions  is  constructed  from  Eq.  (3)  in  which  the  two  row  vectors  of 
atomic  eigenstates  and  are  transformed  to  eigenstates  quantized  in  and 

functions  of  a  co-aligned  2-axis  coordinate  system  in  which  diatomic  calculations  are  con¬ 
veniently  performed,  with  ah  other  functions  unchanged  [28], 

*S,/3)(rd  :  R)  =  {$(1)(i)  0  •  •  •  \id )  0  •  •  •  ^\jd)  0  •  •  •  $(Ar)(n)}o 

=  $(r  :  R)  ■  jl^  0  ■  -  0  •  ■  0  ■  ■  -  IW} 

=  #(r  :  R)  ■T>{a^\Ra/3^.  (19) 

Here,  $£}(zd)  =  &a\i)  ■  D^(Rapy  and  ^\jd)  =  &W(j)  •  D(/3)(^)t  are  the  (a, /3) 
atomic  spectral  states  in  the  diatomic  frame,  D (a\Rap)  and  D ^\Rap)  are  matrices  com¬ 
prising  appropriately  ordered  sets  of  rotation  matrices  for  ah  of  the  multiplet  eigenstates 
of  the  atoms  a  and  /?,  respectively,  the  Euler  angles  Rap  are  given  by  the  polar  angles 
((f) a /3,  Oap,  0)  specifying  the  direction  of  atom  (3  from  a  in  the  laboratory  frame  [16],  and 
the  matrix  D ^a^(Rap)  defined  by  Eq.  (19)  is  expressed  as  an  appropriately  ordered  outer 
product  of  the  two  sets  of  atomic  rotation  matrices  and  unit  matrices  (T  T )  for  the  other  N- 
2  atoms  [8].  The  notation  rd  employed  in  ^^’^\rd  :  R),  and  in  related  functions  defined 
below,  is  meant  to  indicate  a  change  of  angular  coordinates  ((f) d  =  <f>  —  4>ap,  @d  —  9  ~  ®ap) 
only  for  those  electrons  (i,j  — >  id,jd)  assigned  to  the  atoms  (a,  (3)  whose  eigenstates  have 
been  transformed  to  the  diatomic  coordinate  frame,  with  the  other  nt  —  na  —  rip  electron 
coordinates  and  IV  —  2  row  vectors  of  atomic  states  represented  in  the  laboratory  frame. 

Transformation  of  the  diatomic  Hamiltonian  of  Eq.  (16)  to  the  basis  of  Eq.  (19)  provides 
the  representation 

H {a'P\Rap)  =  D {a^\Rap)]  ■  H ^(Rap)  •  B^\Rap),  (20) 

where  the  diatomic  Hamiltonian  matrix 

H ^(Rap)  =  :  R)\H^\id,jd  :  Rap)\^\rd  :  R))  (21) 
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depends  by  construction  only  on  the  scalar  separation  Rap.  That  is,  the  calculations 
implied  by  Eq.  (21)  are  performed  in  the  co-aligned  diatomic  coordinate  systems  in  which 
the  scalar  Rap  specifies  the  separation  of  the  atoms  a  and  (3  along  the  common  2-axis. 


2.2.2  Diatomic  Subspace  Representations 

The  physical  and  unphysical  blocks  of  the  Hamiltonian  of  Eq.  (21)  can  be  isolated  fol¬ 
lowing  a  development  similar  to  that  of  Eqs.  (8)  to  (11).  In  this  diatomic  case,  however, 
transformation  to  a  permutation-symmetry-adapted  representation  is  required  only  in  the 
subspace  of  product  functions  Cg>  &^\jd)  associated  with  the  atoms  (a,  (3)  in 

which  the  diatomic  Hamiltonian  of  Eq.  (14)  acts,  although  it  is  convenient  to  continue 
to  write  out  all  matrices  in  the  dimensions  of  the  full  spectral-product  basis  of  Eq.  (3). 
Following  this  convention,  the  unitary  transformation  [cf.,  Eq.  (8)] 

H (s’0\Ra0)  =  U •  V(sa'n(Ra0) 

=  ( {u‘S’0\Rd,)}m  /{H tP\Rc.p)}pp  0  \ 

(22) 

is  employed  to  isolate  the  physical  and  unphysical  blocks  of  the  diatomic  Hamiltonian  in 
the  co-aligned  coordinate  frame,  where  Ug“’,8')(.RQ,/3)  is  the  unitary  matrix  obtained  from 
the  diagonalization  [cf.,  Eqs.  (9)  and  (10)] 


pj(a>d) 

{S 


s  (Red)'  ■  S<g-I3) (Ral))  ■  =  s\‘ 


A) 


( Ra/3 ) 


(a,/3)( 
d 


o 


n(a.d)  I 

-*-pp 

0 


(23) 


o 


(Rap)}pp 

of  the  Hermitian  matrix  representative  of  the  atomic-pair  antisymmetrizer  in  the  spectral 
product  basis  of  Eq.  (19), 


o(a,/3) 


(Rap)  =  (Q{a^)l,2(^\rd  :  R)|Pia’/3)|#^)(rd  :  R)) 
=  {P^P)^p\rd  :  R)\P^]^\rd  :  R)). 


(24) 


The  asymptotically  wave-function-normalized  antisymmetrizer  employed  in  Eqs.  (23)  and 
(24),  P^a,/3)  =  ((na  +  np)\)-l/2(na\np\)-1/2  J2p=id}jd(-l)SpPP,  operates  only  on  electrons 
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in  the  sets  id  and  jd ,  the  metric  matrix  S ^’^(Rap)  depends  explicitly  on  Rap  through 
the  overlap  of  the  two  sets  of  atomic  eigenstates  on  atoms  a  and  (3,  the  value  of  the 
redundancy  factor  Qla’d)  —  (na  +  n(g)!/(nQ,!nJg!)  of  the  explicitly  antisymmetrized  basis 
p(«,/3)^(«,/3)(rrf  .  follows  from  wave- function  normalization  of  the  antisymmetrizer,  and 
Sq'  j\Raf3  —■ *■  oo)  — >  I  in  the  diatomic  dissociation  limit.  As  in  Eq.  (9),  those  eigenvalues 
in  the  diagonal  matrix  S ^^\Rap)  which  approach  Q are  ordered  in  decreasing  value 
along  the  upper  diagonal  and  those  that  approach  zero  are  placed  in  the  lower  diagonal, 
with  the  corresponding  eigenvectors  of  the  S^,f3\Rap)  matrix  appearing  in  appropriate 
columns  of  the  transformation  matrix  Uga’^(i?ayg)  of  Eq.  (23). 

Of  course,  the  expressions  of  Eqs.  (22)  to  (24)  include  the  matrix  element  re-arrangements 
required  to  enforce  the  ordering  convention  of  Eq.  (3),  signified  by  the  indicated  subscripts, 
resulting  in  matrices  which  disguise  the  blocked  structures  of  the  explicit  expressions  given 
above  [8,12],  That  is,  although  the  eigenvalues  of  the  metric  matrix  of  Eq.  (24)  can  be 
ordered  as  indicated  above  in  matrices  representing  the  diatomic  subspace  of  the  spectral- 
product-  basis,  their  positions  in  the  matrices  of  Eq.  (23)  must  be  in  accordance  with  the 
particular  ordering  convention  adopted  in  Eq.  (3).  Similarly,  the  energy  matrix  of  Eq. 
(22)  will  not  exhibit  the  explicitly  blocked  structure  shown  there  after  the  matrix  elements 
have  been  re-ordered  in  compliance  with  the  ordering  convention  of  Eq.  (3). 

States  of  tot  ally- antisymmetric  and  non-totally-antisymmetric  permutation  symmetry  in 
electron  coordinates  id  and  j,i  defining  the  diatomic  Hamiltonian  matrix  Hga’^  (i?a(g)  of 
Eq.  (22)  take  the  partitioned  form  [cf.,  Eq.  (11)] 

$^>/3)(rd  :  R)  =  ^ \rd  :  R)  ■  Xj£’P\Raf3) 

=  $(r  :  R)  ■  B^\Raf3y  •  \j£’p)(Ral3) 

-  {{*s,P\rd  ■■  R)}p,  (4Q,/3)(^  :  (25) 

where  Eq.  (19)  has  been  employed  in  the  second  line  and  closure  implied  in  the  last  line. 
Finally,  the  dimensions  of  the  physical  subspaces  defined  in  Eqs.  (22)  to  (25)  for  the 
different  atomic  pairs  can  generally  differ  from  one  another  and  from  the  dimension  of  the 
aggregate  physical  subspace  of  Eq.  (11),  as  is  further  clarified  in  Section  2.3  below. 
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2.2.3  Diatomic  Eigenspectra 


The  eigenstates  and  energies  of  the  partitioned  diatomic  Hamiltonian  of  Eq.  (22)  are 
obtained  in  the  usual  way  from  solution  of 


R^\Rap)  •  U ^\Rap)  =  U ^{Rap)  ■  E W\RaP), 


r(a.d) 


(a,/3) , 


where  the  diagonal  eigenvalue  matrix 

and  the  unitary  transformation  matrix 

u  = 


PP 


0 


{E  ^\Rafi)}r 


o 


{ug’«(  R^))m 

q  JTT(a’d) 


0 


(26) 


(27) 


(28) 


{Ug’p\RaP)}uu/  0 

comprising  the  diatomic  eigenvectors  are  block  diagonal  by  construction.  The  associated 
diatomic  eigenstates  take  the  partitioned  form 


^a^{rd  :  R)  =  ^\rd  :  R)  •  U {^\Rap) 

=  ^P\rd  :  R)  ■  U ^(Rap)  ■  V^P\Raf3) 

=  $(r  :  R)  ■  D {a^\kap)]  ■  (Raf3)  •  U {^\Rap) 

^  {{^\rd  :  R)}p,  {^\rd  :  R)}u}0,  (29) 

consequent  of  Eqs.  (25)  to  (28).  Following  the  standard  convention  [17],  the  eigenvalues 
of  the  physical  and  unphysical  blocks  of  Eq.  (27)  are  ordered  in  increasing  value  down 
the  diagonals,  and  the  associated  eigenstates  of  Eq.  (29)  are  correspondingly  ordered  by 
appropriate  arrangement  of  the  columns  of  the  two  blocks  of  the  transformation  matrix 
of  Eq.  (28).  As  in  Eqs.  (22)  to  (24),  the  arrangement  of  elements  in  Eqs.  (27)  to  (29) 
are  brought  in  compliance  with  the  indicated  ordering  convention,  disguising  the  blocked 
natures  of  these  expressions  but  not  eliminating  them. 

Equations  (26)  to  (29)  complete  the  connection  between  the  diatomic  eigenstates  and  the 
aggregate  spectral-product  states,  with  the  three  transformation  matrices  in  Eq.  (29) 
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required  to  relate  the  two  representations.  The  transformation  matrix  Ug*’^  (i?ayg)  of  Eq. 
(23),  obtained  from  the  metric  matrix  of  Eq.  (24),  and  the  matrices  U ^’^{Rap)  and 
E(a>d)  (Rap)  of  Eqs.  (26)  to  (28)  must  be  constructed  individually  for  all  contributing 
atomic  pairs  in  order  to  implement  the  present  development.  These  matrices  depend  only 
on  the  scalar  separation  Rap  and  require  diatomic  calculations  which  are  performed  once 
and  for  all  over  a  suitable  range  of  interatomic  separations,  whereas  the  transformations 
of  Eqs.  (19)  and  (20)  involve  only  the  known  universal  Wigner  D  functions  and  can  be 
evaluated  employing  a  variety  of  available  methods  [29]. 

2.3  Atomic-Pair  Aggregate  Hamiltonian  Matrix 

When  the  atomic  and  diatomic  information  of  the  preceeding  Section  2.2  is  available  from 
previous  calculations,  the  development  of  the  Eqs.  (19)  to  (29)  can  be  inverted  to  give 

H  {s’p\R)  =  Us(JR)t  ■  H  ^\RaP)  •  Us  (R) 

=  ,/3)(i2)f  •  Hjf  ’ ,p\Rccp)  •  Us,t3)(R)  (30) 

for  the  diatomic  Hamiltonian  matrix  of  Eq.  (18),  where  Hg*’^  (i?a/g)  obtained  from  Eqs. 
(26)  to  (28)  is  block  diagonal  by  construction  and  the  transformation  matrix 

U  {s,I3\R)  =  U^l/3)(i?^)t  ■  D  ^\Raf3)  ■  Us  (R)  -  (^^(rd  :  R)\$s(r  :  R))  (31) 

defined  by  Eq.  (30)  is  the  product  of  the  three  indicated  transformation  matrices.  Each  of 
the  latter  matrices  has  the  dimension  of  the  full  spectral-product  basis,  with  the  indicated 
limiting  expression  obtained  employing  the  definitions  of  Eqs.  (11),  (19),  and  (25)  and 
summing  over  the  formally  complete  intermediate  states. 

The  closure  limit  of  Eq.  (31)  demonstrates  that  the  matrix  Ug y’^(R)  has  a  zero  block  in 
its  lower  left-hand  corner,  {Ug“’/;i)(.R)}Up  =  0,  since  {$ga’/3')(r,d  :  R)}u  of  Eq.  (25)  is  non- 
totally  antisymmetric  in  electrons  in  the  sets  id  and  j,i  and  {<I>s(f  :  R)}p  of  Eq.  (11)  is 
totally  antisymmetric  in  all  electrons  in  this  limit.  This  circumstance  is  sufficient  to  insure 
that  the  upper  left-hand  block  {HgQ:/3'1  (R')}pp  of  Eq.  (30)  requires  for  its  construction  only 
the  physical  blocks  of  the  diatomic  energy  and  transformation  matices  of  Eqs.  (27)  and 
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(28),  and  of  the  transformation  matrix  Uga’^(.R)  of  Eq.  (31),  with  contributions  from  the 
unphysical  and  off-diagonal  blocks  of  these  latter  matrices  vanishing  in  the  closure  limit. 

In  light  of  the  foregoing,  the  unphysical  diatomic  solutions  of  Eqs.  (26)  to  (28)  are  not 
required  in  constructing  the  physical  block  of  the  aggregate  Hamiltonian  matrix  of  Eq. 
(17),  which  can  be  written  in  the  three  equivalent  forms 

N  N  /  yy  r)\ 

cy=l  /3=a-\- 1 

N  N  /  j xj 

a=l /3=a+l  ^  ’ 

N  N 

=  2E{{H^>(Ji)}pp+  5E  (32) 

OL  =  1  (3=&-\- 1 

where  the  physical  block  of  the  diatomic  Hamiltonian  matrix  obtained  from  Eqs.  (30)  and 
(31)  is 

{4“’«(R)}pp  =  {U^“ ’«(R)}tp  ■  {H {tn{RMrr  •  {Uga'^(R)}pP.  (33) 

Equations  (32)  and  (33)  make  clear  that  the  dimensions  of  the  physical  subspace  vectors 
(r,i  :  R)}p  and  the  square  physical  blocks  of  the  atomic-pair  Hamiltonian  matrices 
{Hg*’^  (Rap)} pp  can  generally  be  different  for  the  different  atomic  pairs,  and  that  they 
can  also  differ  from  the  dimension  of  the  physical  subspace  { d*  s  ( r  :  R)  }P  of  the  aggre¬ 
gate  representation  of  Eq.  (11).  Accordingly,  the  transformation  matrices  { Ug0’1'^ (i?j }pp 
obtained  from  Eq.  (31)  are  generally  rectangular,  with  dimensions  appropriate  to  the  in¬ 
dicated  matrix  multiplications  of  Eq.  (33),  the  final  matrices  {Hg“’^')(i?)}pp  all  having  the 
dimension  of  the  physical  subspace  of  Eq.  (11). 

The  notation  employed  in  the  second  term  in  the  first  line  of  Eq.  (32)  implies  that  the 
scalar  limit  Rap  — >  oo  is  to  be  taken  in  all  matrices  in  Eqs.  (30),  (31),  and  (33)  except  the 
transformation  matrix  Us(-R).  Accordingly,  this  second  term  is  seen  to  reduce  correctly 
to  the  indicated  sum  of  atomic  Hamiltonian  matrices  required  to  cancel  the  additional 
atomic  terms  included  in  the  first  double  sum  over  atomic  pairs.  Specifically,  employing 
the  limits  E {a^(Rap  ->  oo)  and  U £,P){Rap  — ►  oo)  — ►  I  in  Eq.  (33),  and 
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noting  that  Ug*’^  (i?aig  — »  oo)  — ■»  Ug“’,8')  (oo)  in  Eq.  (31),  gives  U ^,i3\r)  — ■>  Ug*’^  (oo)  • 
n^\Raf3)-us(R).  Finally,  noting  that.  U  g*’^  (OO)  will  only  include  off-diagonal  terms 
which  correspond  to  admixtures  of  degenerate  states  in  the  sum  of  matrices  H  and 
Hhh  by  construction,  and  employing  the  commutation  of  D ^a,P\Rap)  with  the  atomic 
energy  matrices  HrT  and  Hf  d  consequent  of  the  well-known  energy  degeneracy  of  atomic 
multiplet  states  [18],  demonstrates  the  equivalence  of  the  three  forms  of  the  Hamiltonian 
matrix  given  in  Eq.  (32). 

Equation  (32)  and  the  associated  transformations  of  Eqs.  (30),  (31),  and  (33)  provide  a 
formally  exact  expression  for  the  Hamiltonian  matrix  of  a  molecule  or  other  finite  atomic 
aggregate  in  terms  of  atomic  and  diatomic  energies  which  can  be  determined  once  and  for 
all  from  appropriate  separate  calulations.  The  metric  and  transformation  matrices  of  Eqs. 
(9)  and  (10)  are  the  only  factors  in  the  development  that  refer  explicitly  to  the  particular 
spatial  arrangement  of  the  interacting  atomic  constituents.  These  provide  the  modifica¬ 
tions  of  the  bare  atomic  energies  and  diatomic  interaction  energies  due  to  the  inter-atomic 
aggregate  electron  antisymmetry  required  upon  their  incorporation  in  the  larger  aggregate, 
and  must  be  determined  separately  for  each  individual  atomic  arrangement  considered. 
Since  aggregate  electron  antisymmetry  is  expressed  in  the  present  development  through 
the  Us(-R)  matrix,  rather  than  as  an  ad  hoc  global  constraint,  accurate  evaluations  of 
the  interaction  Hamiltonian  of  Eq.  (33)  may  be  obtained  in  some  cases  by  including  only 
small  numbers  of  atoms  in  near  proximity  to  a  given  interacting  pair  in  construction  of 
the  required  metric  and  transformation  matrices.  More  generally,  essentially  complete 
evaluations  of  the  aggregate  metric  matrix  can  be  performed  from  the  one-electron  or¬ 
bital  overlap  integrals  in  which  it  can  ultimately  be  expressed,  and  the  associated  unitary 
transformation  matrix  obtained  from  Eq.  (9). 

The  atomic  and  diatomic  calculations  required  in  constructing  Eqs.  (32)  and  (33)  are 
performed  variationally,  insuring  that  the  resulting  atomic  energies  and  potential  energy 
curves  provide  upper  bounds  to  the  corresponding  energies  obtained  from  exact  solutions 
of  the  associated  atomic  and  diatomic  Schrodinger  equations.  These  calculations  fix  the 
dimensions  of  the  physical  vectors  {$ga,/3') (rd  :  R)}P  appearing  in  Eq.  (31)  and  of  the 
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physical  blocks  {Hga,,S')(i?Q,(g)}pp  of  the  diatomic  interaction  Hamiltonian  matrices  appear¬ 
ing  in  Eq.  (33),  but  not  necessarily  that  of  the  dimension  of  the  aggregate  physical  vector 
{$s(r-  :  R)}P  in  Eq.  (31).  Accordingly,  as  the  dimension  of  the  latter  representation  is 
increased  for  fixed  dimensions  of  the  diatomic  subspaces,  bounded  (Hylleraas-Undheim) 
convergence  is  obtained  to  the  aggregate  energies  defined  by  the  chosen  diatomic  repre¬ 
sentations.  Of  course,  this  does  not  insure  that  the  aggregate  energies  obtained  in  this 
manner  will  provide  upper  bounds  to  those  energies  obtained  from  exact  solutions  of  the 
Schrodinger  equation  of  Eq.  (1).  Nevertheless,  it  can  be  anticipated  that  reliable  results  can 
be  obtained  from  the  development  for  specific  aggregate  states  when  sufficiently  converged 
representations  are  achieved  for  the  corresponding  bare  diatomic  states,  issues  which  must 
be  addressed  in  the  course  of  computational  implementations  and  specific  applications. 

3.0  Illustrative  Computational  Implementation 

A  particular  computational  implementation  of  the  foregoing  development  is  provided  in 
Section  3.1  employing  valence-bond  methods  in  constructing  atomic  and  diatomic  informa¬ 
tion,  explicit  closure  limits  of  the  aggregate  spectral-product  representation  in  completing 
the  evaluation  of  the  exact-pair  expressions  of  Eqs.  (32)  and  (33)  are  described  in  Section 

3.2,  and  an  illustrative  computational  application  to  the  H3  molecule  is  reported  in  Section 

3.3. 

3.1  Calculations  of  Atomic  and  Diatomic  Information 

The  diatomic  information  required  in  evaluation  of  the  matrices  of  Eqs.  (30)  to  (33),  al¬ 
though  formally  defined  in  the  spectral-product  representation,  can  be  obtained  employing 
explicitly  antisymmetrized  valence-bond  methods  [5].  This  is  conveniently  accomplished 
in  the  present  implementation  by  construction  of  an  explicitly  antisymmetric  form  of  the 
permutation-symmetry-adapted  diatomic  subspace  :  Rap)  of  Eq.  (25).  Specifi¬ 

cally  required  for  this  purpose,  in  addition  to  the  atomic  eigenstates  <&(“)(*)  and  energy 
matrices  are  the  diatomic  matrices  S U ^’^(Rap)  and  the  diatomic  eigen- 

spectra  {E^,/3)  (Ra/3)}pp:  {U h ’^')(f?a/3)}pp  °f  Sections  2.2.2  and  2.2.3  obtained  from  the 
Hamiltonian  matrix  {Hga’^\f?a/g)}pp,  all  of  which  information  is  provided  by  suitable  pro- 
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cessing  of  the  results  of  explicitly  antisymmetric  valence-bond  calculations  incorporated 
into  the  spectral-product  development. 

Sets  of  orthonormal  antisymmetric  atomic  configurational  state  functions  4T“)(7),  cho¬ 
sen  in  the  form  of  Slater  determinants  of  canonical  Hartree-Fock  orbitals,  are  employed 
in  construction  of  the  required  row  vectors  of  atomic  spectral  eigenfunctions  $ ia^ (i)  — 
3>(a'(i)-lVTa).  Here,  the  matrix  1VF“)  is  obtained  from  diagonalization  of  the  atomic  Hamil¬ 
tonian  matrix  (4y")(*)|lT“)(i)|<I>^a)(*))  constructed  in  the  orthonormal  basis  (i), 
which  calculations  also  provide  the  atomic  energy  matrices  of  Eq.  (7).  The  ag¬ 

gregate  atomic  spectral-product  basis  of  Eq.  (3)  is  obtained  in  this  way  in  the  form 
T*(r  :  R)  —  $(r  :  R )  •  M,  where  [cf. ,  Eq.  (3)] 

4>(r  :  R)  =  |$(1)(i)  ®  4>(2)(2)  ®  •  •  •  $(A°(n)}o  (34) 

is  the  spectral  product  of  atomic  configurational  state  functions  and  the  matrix  M  is  given 
by  the  outer  matrix  product 

M=  {m^  <g>M(2)  (35) 

Similarly,  the  (a,  /3)  metric  matrix  of  Eq.  (24)  is  obtained  in  the  form 

S ■  S^>/3)(i?a/3)  •  (36) 

where  is  the  outer  matrix  product  [cf.,  Eq.  (35)] 

M(a>/3)  =  |j(l)  0  . . .  m(«)  0  . . .  m(/3)  ®  •  I(Ar) }  .  (37) 

The  metric  matrix  S ^’^(Rap))  —  (Q('a’^)1/2(&i),l3\rd  ■  -R)|-Pjja,,S'>  |3>d '^\rd  :  R))  is 
obtained  from  valence-bond  calculations  of  the  states  [cf.,  Eq.  (19)] 

^P\rd  :  R)  =  |$(1)(i)  ®  •  •  •  ®{u\id)  ®  •  •  •  ®{v\jd)  0  ■  ■  •  $(A°(n)}0  (38) 

employing  Eq.  (34)  and  the  indicated  atomic  spectral  substitutions.  The  information  of 
Eqs.  (34)  to  (38),  when  employed  in  the  development  of  Eqs.  (22)  to  (24)  providing  the 
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unitary  transformation  matrix  Ug*’^  (i?aig)  and  eigenvalues  S ^’^\Rap),  gives  the  physical 
block  of  the  diatomic  Hamiltonian  matrix  of  Eq.  (22)  in  the  familiar  form  [11,12] 


t'-'S 


Here,  ^\Rap)  =  •  Hj ^,l3\Rap)  •  M*^510)  expresses  the  diatomic  Hamiltonian 

matrix  of  Eq.  (21)  in  terms  of  the  corresponding  matrix  Hp  (Rap)  constructed  in  the 
lion-orthogonal  valence-bond  basis  ’^\rd  ■  R),  with  ^^’^\rd  :  R)  given  by  Eq. 

(38).  Although  the  foregoing  expressions  are  presented  in  the  dimension  of  the  aggregate 
spectral-product  basis  as  a  notational  convenience,  of  course,  only  atomic  and  diatomic 
representations  are  required  in  the  calculations  performed. 


Finally,  it  is  seen  from  the  development  of  Eqs.  (34)  to  (39)  that  the  physical  diatomic 
subspace  of  Eq.  (25)  is  obtained  from  the  present  implementation  in  the  form 


{^\rd  :  R)}p  =  :  R )  •  U ^ (Raf3)}P 

=  {P{A^^\rd  ■■  R)  •  U  •  {S  ^(Rap)}^2 

-  {P{Aa^^\rd  :  R)  ■  U {sa^\Ra/3)}P  (Q{a,p)y1/2,  (40) 

where  the  '^{Rap)  and  Ug ?^\Rap)  matrices  and  the  basis  states  ^^^\rd  :  R)  = 
^\rd  '■  R)  •  are  obtained  from  the  foregoing  expressions.  Equations  (34)  to 

(40)  are  evidently  sufficient  to  obtain  the  physical  blocks  of  the  diatomic  energy  and 
transformation  matrices  by  the  diagonalization  of  Eqs.  (27)  and  (28),  which  completes 
determination  of  the  specifically  atomic  and  diatomic  information  required  in  constructing 
the  atomic-pair  interaction  energy  matrix  of  Eq.  (33) . 

The  role  of  the  Ug“’^  (i?ayg)  matrix  of  Eq.  (23)  in  the  first  line  of  Eq.  (40)  is  to  iso¬ 
late  the  totally  antisymmetric  subspace  of  the  (a,  (5)  spectral-product  representation  in 
accordance  with  Eqs.  (22)  to  (24)  [11,12],  whereas  in  the  second  and  third  lines  the 
transformation  provides  a  linearly-independent  orthonormal  subspace  of  the  Q0v-T-fbld 
redundant,  linear ly-dependent,  explicitly  antisymmetrized  ( a ,  j3)  spectral-product  basis 
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pCa,/3)^(a,/3)(^  .  constructed  following  well-known  canonical  orthogonalization  proce¬ 
dures  [30].  Correspondingly,  the  spectral-product  expression  for  (r,/  :  R)}p  in  the 

first  line  in  Eq.  (40)  converges  to  the  indicated  antisymmetric  expression  for  diatomic 
separations  Rap  at  which  the  spectral  closure  limit  is  attained  [11,12],  As  a  consequence, 
the  equivalence  between  the  spectral-product  form  and  the  explicitly  antisymmetrized  ex¬ 
pressions  of  Eq.  (39)  can  fail  in  the  dissociation  limit,  when  S^’^\Rap  —■ >  oo )  — >  I  rather 
than  — >  Q&’ftl,  with  the  spectral-product  representation  of  t^\rd  '■  R)}p  dissociating 
to  a  product  form  and  the  explicitly  antisymmetric  expressions  in  the  second  and  third 
lines  of  Eq.  (40)  remaining  antisymmetric  in  the  dissociation  limit. 

3.2  Aggregate  Spectral-Product  Closure  Limits 

Completion  of  construction  of  the  Hamiltonian  matrix  of  Eqs.  (32)  and  (33)  requires  the 
aggregate  physical  subspace  {ffisA'  :  R)}p  of  Eq.  (11)  and  the  transformation  matrix 
of  Eq.  (31).  Following  the  development  of  Eqs.  (34)  to  (40),  the  aggregate  subspace  is 
obtained  in  the  form 

{$s (r  :  R)}P  =  {$(r  :  R)  •  U S(R)}P 

=  {PA*(r  :  R)  ■  Us(fi)}p  ■  {Sd(fl)};y2 

-  {Pa’Ht  :  III  ■  US(B)}P  Q-1^,  (41) 

which  expressions  have  been  shown  to  converge  as  indicated  for  all  atomic  arrangements 
R  at  which  the  limit  of  spectral  closure  is  achieved  [11,12],  Here,  the  aggregate  basis 
#(r  :  R)  =  &(r  :  R)  M  is  obtained  from  Eqs.  (34)  and  (35)  and  the  transformation 
matrix  Ug(-R)  from  Eqs.  (9)  and  (10)  employing  the  metric  matrix  constructed  in  the 
valence-bond  representation,  S (R)  =  (51/2M^’  ■  (4>(r  :  R)\pA\&(r  :  R))  ■  M.  Although 
the  second  and  third  lines  of  Eq.  (41)  as  written  formally  imply  global  enforcement  of 
antisymmetry  for  all  aggregate  electrons,  when  the  range  of  interactions  in  the  Hamiltonian 
matrix  of  Eqs.  (32)  and  (33)  are  sufficiently  limited  only  small  numbers  of  atoms  in  close 
spatial  proximity  to  a  given  interacting  pair  will  contribute  significally  to  the  metric  matrix 
Us(-R)  required  for  evaluation  of  Eq.  (33),  potentially  simplifying  the  calculation.  More 
generally,  as  indicated  above,  complete  evaluation  of  the  transformation  matrix  Us(-R)  is 
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accomplished  from  the  one-electron  overlap  integrals  which  are  required  in  the  construction 
of  the  corresponding  S-matrix. 

The  Us(-R)  matrix  in  the  first  line  of  Eq.  (41)  isolates  the  totally  antisymmetric  subspace 
of  the  spectral-product  representation,  in  accordance  with  Eqs.  (8)  to  (11),  whereas  in 
the  second  and  third  lines  it  provides  a  linearly-independent  orthonormal  subspace  of  the 
<3 -fold  redundant,  linear ly-dependent,  explicitly  antisymmetrized  spectral-product  basis 
[11,12],  As  in  Eq.  (40),  the  equivalence  of  the  spectral-product  expression  in  Eq.  (41) 
with  two  explicitly  antisymmetrized  expressions  given  there  can  fail  in  the  dissociation 
limit  when  S (R  — >  oo)  — >  I  rather  than  — >  Q I.  Correspondingly,  as  indicated  in  the 
discussion  of  Eq.  (12)  above,  the  spectral-product  representation  of  {$s(^  :  R)}p  in  this 
case  dissociates  to  a  product  form  [<I?(r*  :  oo)  •  Us(oo)],  whereas  the  expressions  in  the 
second  and  third  lines  of  Eq.  (41)  remain  explicitly  antisymmetric  in  the  dissociation  limit 
in  the  form  Pa^{v  :  oo)  •  Us(oo). 

Finally,  employing  Eqs.  (40)  and  (41)  in  evaluating  the  physical  block  of  the  transformation 
matrix  Ug “’^(J?)  of  Eq.  (31)  gives 

{l4“’«(K)}„„  =  {(44“-»(rd  :  H)|#s(r  :  R)))m 

=  {S (P\R^)IQ{O-0XP  •  ■  D<“-»(fi„y  ■  Us(fi)}„  •  {S dR)/Q}X? 

-  { •  IV'V/i,.,,)  •  Us (/()}, „„  (42) 

where  the  last  line  indicates  the  closure  limit.  The  latter  expression  is  evidently  identical  in 
appearance  to  the  physical  block  of  Eq.  (31)  constructed  entirely  in  the  spectral-product 
representation,  in  spite  of  the  fact  that  Eq.  (42)  employs  the  explicitly  antisymmetric 
diatomic  and  aggregate  states  of  Eqs.  (39)  and  (41).  This  equivalence  is  in  accord  with 
the  recognition  that  any  distinction  between  the  physical  subspace  of  the  atomic  spectral- 
product  representation  and  the  linearly-independent  subspace  of  its  explicitly  antisym¬ 
metrized  form  must  become  inconsequential  in  the  closure  limit  [11,12], 

Although  it  can  be  expected  that  the  explicitly  antisymmetric  expressions  in  Eqs.  (40) 
and  (41)  will  be  more  rapidly  convergent  with  dimension  of  the  basis  than  will  be  the 
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expansions  in  the  product  basis  itself,  some  care  must  be  exercised  in  using  the  former 
expressions.  Specifically,  the  aforementioned  failure  of  equivalence  between  the  spectral- 
product-  and  explicitly  antisymmetric  representations  in  Eqs.  (40)  and  (41)  at  arbitrarily 
large  Rap  and  R  results  in  the  second  line  of  Eq.  (42)  having  the  dissociation  limit 
{U ^\R  -  oo)}pp  -  (Q^/Q)1/ 2{u£*,/J)(oo)  •  Us(oo)  •  D ^\Ra/3)}pp,  which  in¬ 
cludes  the  apparently  spurious  factor  (Q(a’^ /Q)1/2.  This  factor  would  seem  to  mod¬ 
ify  the  long-range  behavior  of  the  matrix  {HgV’^(fl)}J)))  of  Eq.  (33),  which  has  the 
apparent  limit  {H {°>f3)(R  oo)}pp  (Q^/Q)  Us(oo)t  •  {H^  +  H^)}  •  Us(oo), 

causing  Eq.  (32)  to  violate  the  dissociation  limit  of  Eq.  (12).  As  a  practical  mat¬ 
ter,  however,  the  expressions  implied  by  the  last  lines  of  Eqs.  (40)  to  (42)  can  always 
be  satisfied  in  large  but  finite  spatial  regions.  Specifically,  the  spectral  closure  limits 
<  R^)}„  -  Q<“'« Ipp  can  be  satisfied  for  each  atomic  pair  in  the  indicated 
intervals  in  suitably  chosen  denumerable  L2  diatomic  representations,  where  the  radii 
have  sufficiently  large  but  hnite  values  such  that  the  diatomic  energy  matrices  in  Eq.  (33) 
take  the  limiting  forms  {^a^\Rap  — >  R-[^  ) }pp  — >  H(a  ^  +  H*^.  The  spectral  closure 
limit  {S d(R  <  Roo)}PP  — ?  Q1  is  correspondingly  attained  in  the  domain  R  <  R^o  employ¬ 
ing  the  aggregate  basis  states  defined  by  at  least  the  individual  diatomic  representations 
employed,  where  R0 0  is  the  union  of  the  ^  values  for  all  atomic  pairs.  Under  these 
conditions,  the  closure  limits  implied  in  the  last  lines  of  Eqs.  (40)  to  (42)  are  satisfied  and 
the  aggregate  Hamiltonian  matrix  {Hs(-Roo)}pP  of  Eq.  (32)  can  be  correctly  regarded  as 
the  sum  of  non-interacting  atomic  Hamiltonian  matrices  of  Eq.  (12). 

3.3  Illustrative  Computational  Application 

The  doublet  and  quartet  potential  energy  surfaces  of  the  H3  molecule  are  studied  in  a 
minimal  (Is3)  basis  [31]  as  an  illustrative  application  of  the  particular  implementation  of 
the  exact-pair  formalism  given  in  the  foregoing  Sections  3.1  and  3.2.  Of  course,  closure 
limits  are  not  achieved  in  this  model  representation,  which  provides  only  an  illustration  in 
which  it  is  possible  to  write  out  in  detail  the  various  expressions  devised,  to  demonstrate 
that  sensible  results  can  be  obtained  from  the  development  employing  even  the  smallest 
possible  basis  sets,  and  to  aid  comparisons  with  related  semi-empirical  approaches  [13,14], 

This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited.  Page  27 


Restricting  attention  to  Ms  —  +1/2  states  of  the  H3  molecule,  the  spectral-product  basis 
of  Eq.  (3)  in  this  case  takes  the  three-term  form 

$(r  :  R)  =  (a(l)a(2)/?(3),a(l)^(2)a(3),/?(l)a(2)a(3)}lsa(l)ls6(2)lsc(3),  (43) 


where  r  =  (r  1, 7*2,  7*3),  R  =  (i2a  =  0,  Rbl  Rc),  a(i)  and  j3{i)  are  the  Pauli  spin  functions, 
and  lsa(i)  =  lSa(|7*i  —  -Ra|)  refers  to  the  Is  orbital  centered  on  atom  a.  The  metric  matrix 
of  Eq.  (10)  in  this  case  is 


S  (R) 


1  -  s2 

1  bab 

SabSacSbc  Sbc 
Sab^acSbc  ®ac 


^ab^ac^bc  $bc 


SabSacSbc  ^ab 


SabS 

SabS 


acsbc  sac 
2 

acSbc  ~  Sab 


(44) 


where  sap  —  s(Rap)  =  (lsa(i)|ls/3(i))  <  1  is  the  one-electron  Is  orbital  overlap  integral, 
and  the  transformation  matrix  Us(-R)  and  eigenvalue  matrix  S d(R)  are  obtained  from 
diagonalization  of  the  metric  matrix  of  Eq.  (44)  following  the  development  of  Eqs.  (9)  and 
(10).  Specifically,  the  eigenvalues  of  the  S(J2)  matrix  are  found  to  be 


{Sd(-R)}l,2  —  1  sabsbcsca  i  ^y^{(Sab  sac)~  T  (sab  Sbc )  +  (sac  Sbc )  }  (45a) 


{Sd(fi)}3  =  1  -si,- 

& ac  $bc  T  2 Sab^bc^aci 

(45  6) 

and  the  corrersponding  Us  (R)  matrix  is 

(  76COS0(jR)  “  75 sinQ{R) 

-^sinQ(R)  +  -^cosQ(R) 

Us(-R)  = 

-^cosQ(R)  +  -^sinQ(R) 

-±=sinG{R)  -  -±=cosG{R) 

i 

V3 

,  (46a) 

\  ~^cos®(R) 

~7E sin®(R. ) 

73  / 

where 

Q(R)  =  — tan _1  j 

V3(s2ab  -  s2ac)  \ 

(466) 

(4b  s2bc  +  s2ac)r 

The  two  eigenvalues  given  by  Eq.  (45a)  and  the  first  two  columns  of  the  matrix  of  Eq. 
(46a)  correspond  to  doublet  states,  whereas  Eq.  (45b)  and  the  last  column  of  Eq.  (46a) 
correspond  to  a  quartet  state.  That  is,  the  states  $s(^  :  R)  of  Eqs.  (11)  and  (41)  in  this 
example  are  ,5=1/2,  3/2;  M$=+ 1/2  spin  multiplets,  symmetries  which  are  maintained  in 
the  three-atom  dissociation  limit,  in  which  case  cosQ(R  — »  oo)  — >  1,  sinQ(R  — >  oo)  — > 
0,  and  Us  (R  — >  oo)  — >  Us(oo)  7  I.  They  are  linearly-independent  combinations  of 
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the  antisymmetrized  forms  of  the  three-term  basis  of  Eq.  (43)  when  the  corresponding 
eigenvalues  of  Eqs.  (45)  do  not  vanish. 

The  diatomic  metric  matrices  ’^(Rap)  and  the  corrersponding  transformation  matrices 
Uga,/3')  (Rap)  of  Eqs.  (23)  and  (24)  are  obtained  in  the  forms 


S  {vb\Rab)  = 


a(a?c)  (  td  \  _ 

y-^ac)  — 


Sg’c)(i?bc)  = 


0 

0 


1 
0 

-si 


—s 


0 


0 

1 

o 

—s 


ab 


1  -st 


be 


—s 

l 

o 


be 


-K 

o 

l 


0 

0 

1  -  s2 
1  bbc  . 


u 


(a,b) 


( Rab )  = 


U^c)(i?ac)  = 


/  7  \  V  ^  V  A 

u?'c)(jfc)=  -4,  o  4, 


0 

1 

0  > 

1 

\/2 

0 

1 

\/2 

J  ,  (47a) 

1 

0 

JL  j 

V2 

■  1 

0 

\/2 

V2 

\ 

0 

1 

0 

,  (476) 

1 

0 

1 

. 

) 

1 

0 

J_\ 

V2 

1 

\/2 

0 

75 

,  (47c) 

0 

1 

0  / 

with  the  eigenvalues  of  the  three  diatomic  S-matrices  all  taking  the  identical  form 


s {tP)(Rap)  =  u ■  S^\Raf3)  ■  U 

/I  +  4/J  0  0  \ 

=  0  1  -  sip  0  .  (47c/) 

V  o  o  1  -sip) 

Construction  of  the  states  3>g ?,f3\rd  :  R )  of  Eqs.  (25)  and  (40)  employing  the  Ug y,f3\Rap) 
matrices  of  Eqs.  (47a)  to  (47c)  shows  them  to  be  products  of  diatomic  states  (ab,  ac,  or  be) 
°fl£+  or  3E+  symmetry  with  a  Is  atomic  spin  orbital  (c,  b ,  or  a)  of  appropriate  electron 
spin.  That  is,  the  transformations  of  Eqs.  (25)  and  (40)  provide  the  correct  diatomic 
eigenstates  for  the  three  atomic  pairs,  so  that  the  diatomic  Hamiltonian  matrices  of  Eqs. 
(22)  and  (39)  are  diagonal  and  U (Rap)  =  I  for  all  three  pairs  in  this  example.  The 
diagonal  diatomic  energy  matrix  E (Rap)  of  Eq.  (27)  is  the  same  for  all  three  diatomic 
pairs  in  this  case  and  takes  the  form 

/ 1E(Ra/3 )  0  0  \ 

E  {a^(Rai3)=\  0  3E(Raf3)  0  ,  (48) 

\  0  0  3E(Ra/3) ) 


where  1,sE(Rap)  refer  to  the  singlet  and  triplet  Heitler-London  H2  energy  curves. 
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The  aggregate  atomic-pair  Hamiltonian  matrix  of  Eq.  (17)  in  this  case  is 

HS(.R)  =  H (s’b)(R)  +  H {s,c)(R)  +  Hg ,c)(i2)  -  (H^o)(fl)  +  H §\r)  +  H ^{R)} 

=  H  {s,b)(R)  +  H  {s,c)(R)  +  H  {s,c\R)  -  3ElsI 

=  3ElsI  +  {H ^(R)  -  2ElsI}  +  (Hga,c)(i2)  -  2ElsI}  +  {H§’c\r)  -  2ElsI},  (49) 

where  E is  is  the  Is  orbital  energy  and  the  diatomic  pair  terms  of  Eqs.  (30),  (31),  and 
(33), 

Hg°,/3)(JJ)  =  Us(fl)'  ■  ■  Us(-R),  (50) 

are  constructed  from  Eqs.  (39)  to  (42).  Note  that  D ^a,^(Rap)  =  I  in  Eqs.  (49)  and  (50), 
since  only  Is  orbitals  are  employed,  ’^(-Ra/j)  =  I  since  the  correct  Heitler-London  1  E,j 
and  3E+  wave  functions  are  obtained  from  the  permutation-symmetry-adapted  diatomic 
states  of  Eqs.  (25)  and  (40)  in  this  case,  and  the  last  two  lines  in  Eq.  (49)  follow  from  the 
commutation  of  the  three  atomic  Hamiltonian  matrices  (H,fd  =  EisT)  with  the  unitary 
transformation  matrix  Ug(-R).  Finally,  the  aggregate  Hamiltonian  matrix  of  Eq.  (49)  can 
be  written  in  a  form  that  emphasizes  the  separate  contributions  from  atomic  energies  and 
atomic-pair  interaction  potentials, 

Hs(-R)  =  3ElsI  +  Us (R)j  ■  V(ii)  •  Us(fl),  (51) 


Here,  V^v)  =  [3U(RQ/3)  +1V(Ra/3)\/2  is  the  average  and  V^ff)  =  [3V(Raf3)  -  V(RQ/3)]/2 
one-half  the  difference  of  the  triplet  and  singlet  Heitler-London  H2  interaction  energy 
curves,  both  of  which  vanish  in  the  limit  Rap  — >  00. 


Carrying  out  the  indicated  matrix  multiplications  in  Eq.  (51)  shows  that  (Hs(-R)}  13  = 
{Hs(-R)}23  =  0,  as  expected,  and  that  {Hs(-R)}i2  is  generally  small,  vanishing  identically 
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iii  high-symmetry  atomic  arrangements  \Civ-,  Dooh,  D^h}-  In  these  cases,  0(i?)  — >  0  and 
the  diagonal  Hamiltonian  matrix  elements  of  Eqs.  (51)  and  (52)  take  the  simple  forms 

{Hs(R)hi  =  3 Ela  +  i{3  Vo6  +3Vab}  +  i{3  'Vac  +3Vac}  +  Cvbc  +3Vbc}  (53a) 

{Hs(-R)}22  =  3Eis  +  ^{3  3Vab  +1Vab}  +  ^{3  3 Vac  +1KC}  +  ^{3Vbc  +  4  V6c}  (536) 

{Hs(jR)}33  =  3  Els  +  |{3l/a6  +3 14C  +3Hc},  (53c) 

although  more  generally  these  elements  depend  explicitly  on  the  Q(R)  value  of  Eq.  (46b). 
The  particular  energy  values  given  by  Eqs.  (53)  are  helpful  in  examining  the  limiting 
behaviors  of  the  surfaces  in  high  symmetry  and  prove  useful  in  comparisons  with  related 
theoretical  approaches  made  in  the  following  Section  4. 

A  series  of  numerical  evaluations  of  the  foregoing  expressions  are  performed  in  a  12-term 
Gaussian-orbital  representation  of  the  Is  atomic  hydrogen  orbital.  This  basis  provides 
an  orbital  energy  (-0.4999  cm)  in  good  agreement  with  the  correct  value,  and  a  molecular 
hydrogen  ground-state  equilibrium  interatomic  separation  (Re  =  1.642  a o)  and  total  energy 
(-1.1157  au)  in  good  agreement  with  the  accepted  Heitler-London  values  [17].  Conventional 
Gaussian-based  computational  methods  are  employed  in  evaluations  of  all  information 
required  in  this  particular  implementation  of  the  atomic-pair  development  for  H3  [5,6]. 

In  Fig.  1  are  shown  contour  plots  of  the  two  doublet  eigenvalues  of  the  S(_R)  ma¬ 
trix  of  Eqs.  (44)  and  (45)  for  H3  obtained  in  the  Gaussian  Is3  representation  for  noil- 
symmetric  collinear  atomic  arrangements.  The  { S,/(i?.j }  j  eigenvalue  surface  of  Fig.  1(a) 
has  a  maximum  near  the  saddle  point  ( Rs —  x.xx  a o)  of  the  surface,  decreases  rapidly  to 
zero  in  the  united-atom  limit  Rab  =  Rac  — >  0,  approaches  1  in  the  three-body  break-up 
limit  Rab  =  Rac  — >  oo,  and  goes  to  2  in  the  diatomic  limits  Rab  — >  0,  Rac  — >  oo  and 
Rab  — >  oo,Rac  — >  0.  The  vanishing  of  the  eigenvalue  of  Fig.  1(a)  as  Rab  and  Rac  — >  0, 
corresponding  to  formation  of  a  non-Pauli  Is3  atomic  Li  configuration,  indicates  that  the 
minimal  Is3  basis  does  not  support  a  physically  significant  H3  state  near  the  tri atomic 
united-atom  limit.  The  smaller  doublet  eigenvalue  { ( -R) }2  surface  shown  in  Fig.  1(b) 
is  evidently  less  structured  than  that  of  Fig.  1  (a)  and  vanishes  in  both  the  diatomic  and 
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tri atomic  united-atom  limits,  as  does  the  quartet  eigenvalue  (not  shown),  which  limits  can 
be  inferred  from  Eqs.  (44)  and  (45).  The  S-matrix  eigenvalue  surfaces  evidently  identify 
atomic  configurations  in  which  improper  states  arise,  provide  information  complementary 
to  the  associated  energy  surfaces,  and  can  anticipate  the  presence  of  crossings  in  the  cor¬ 
responding  potential-energy  surfaces  on  basis  of  their  degeneracies,  as  demonstrated  in 
further  detail  below. 

Figure  2  provides  constant-energy  contour  plots  of  the  matrix  elements  {Hs(-R)}n  and 
{Hs(-R)}22  obtained  from  Eqs.  (51)  and  (52)  for  the  lion-symmetric  collinear  atomic  ar¬ 
rangements  of  Fig.  1.  These  well-known  energy  surfaces  are  presented  here  to  illustrate  the 
extent  of  agreement  obtained  with  corresponding  conventional  valence-bond  calculations 
in  a  Is3  basis,  to  demonstrate  their  general  accordance  with  accurate  quantum-chemistry 
calculations  [32-35],  and  for  comparison  with  the  less  familiar  S-matrix  eigensurfaces  of 
Fig.  1.  Specifically,  the  calculated  saddle-point-  energy  (-1.5735  au)  of  Fig.  2(a)  compares 
well  with  that  obtained  from  valence-bond  calculations  (-1.5834  au)  in  the  same  basis,  plac¬ 
ing  the  saddle  point  approximately  0.0735  au  (2.00  eV)  below  the  three-body  (H+H-f-H) 
breakup  energy  of -1.500  au  (-40.82  eV),  and  providing  approximately  50%  of  the  accepted 
accurate  saddle-point  binding  energy  of  0.1570  au  (4.273  eV)  relative  to  the  three-body 
breakup  limit  at  a  geometry  (Rs  =  x.xx  ao)  in  general  accord  with  the  accepted  posi¬ 
tion  ( Rs  =  1.757  ao)  [32-35].  The  calculated  saddle-point-  energy  relative  to  the  two-body 
break-up  energy  (H  +  H2)  provides  a  (1.14x  eV)  barrier  to  the  exchange  reaction  which 
is  in  accord  with  the  valence-bond  calculations  (0.876x  eV)  but  is  approximately  a  factor 
of  two  to  three  times  larger  than  the  accepted  value  of  0.4166  eV  [32-35],  emphasizing 
the  well-known  sensitivity  of  this  energy  difference  to  the  computational  approximation 
employed. 

The  structure  of  the  H3  ground-state  energy  surface  of  Fig.  2(a)  is  evidently  complemen¬ 
tary  t-o  that  of  the  associated  S-matrix  eigenvalue  surface  of  Fig.  1(a),  with  the  topography 
of  the  S-matrix  eigensurface  largely  anti-correlated  with  that  of  the  corresponding  energy 
surface,  the  larger  energies  associated  with  smaller  S-matrix  eigenvalues.  Similarly,  the 
energy  surface  for  the  excited  doublet  state  shown  in  Fig.  2(b)  is  seen  to  be  monotonically 
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increasing  from  the  three-atom  break-up  limit  as  the  diatomic  and  tri atomic  united-atom 
limits  are  approached,  and  consequently  anti-correlated  with  the  monotonically  decreasing 
form  of  the  associated  S-matrix  eigenvalue  surface  of  Fig.  1(b).  Of  course,  such  behaviors 
can  be  expected  more  generally  only  for  small  basis  set  evaluations,  or  of  the  lower-lying 
S-matrix  eigensur faces,  in  view  of  the  ultimate  convergence  of  the  highest-lying  eigensur- 
faces  to  the  redundancy  factor  Q  of  Eq.  (9)  in  a  manner  largely  independent  of  atomic 
arrangement.  The  quartet-state  energy  surface  is  found  to  be  repulsive  with  decreasing 
atomic  separations  for  all  non-symmetric  collinear  atomic  arrangements,  in  accordance 
with  the  approximation  of  Eq.  (53c),  and  is  not  depicted  here. 

In  Fig.  3  are  shown  comparisons  of  the  first  two  doublet  energies  in  H3  for  symmetric 
collinear  ( Rab  —  Rac )  atomic  arrangements  obtained  from  Eqs.  (53a)  and  (53b)  with  cor¬ 
responding  valence-bond  calculations.  Evidently,  the  present  results  for  the  ground-state 
energy  surface  are  in  good  agreement  with  the  corresponding  Is3  valence-bond  calcula¬ 
tions  which  include  only  covalent  structures  (o),  and  are  also  in  accord  with  valence-bond 
results  which  include  all  possible  ionic  configurations  (•)  in  the  Is3  basis.  The  first  excited 
doublet  state  in  H3  is  seen  to  be  repulsive,  in  good  agreement  with  the  valence-bond  calcu¬ 
lations  and  in  qualitative  accord  with  previously  reported  accurate  calculations,  although 
this  excited  doublet  is  known  to  undergo  an  avoided  crossing  with  a  higher-lying  state 
which  is  not  included  in  the  present  development  [34],  Although  the  results  of  Fig.  3  for 
the  ground-state  energy  surface  appear  to  bound  the  valence-bond  results  from  above,  this 
is  only  a  happenstance  in  view  of  the  small  basis  set  employed.  Indeed,  similar  apparent 
boundedness  is  not  found  for  the  excited  doublet  state  in  Fig.  3,  although  both  results 
are  clearly  in  satisfactory  agreement  with  both  forms  of  the  valence-bond  calculations, 
the  ionic  configurations  making  only  small  contributions  to  the  ground-state  surface  but 
somewhat  more  significant  contributions  to  the  excited-state  surface. 

Figure  4  depicts  contour  plots  of  the  two  doublet  S-matrix  eigenvalues  in  H3  of  Eq.  (45a) 
for  C2V  arrangements  in  which  the  Rt,i,  —  Rac  separations  are  co-varied  and  the  apex 
angle  6  ranges  from  30  to  180  degrees.  A  seam  of  intersection  of  these  two  surfaces  is 
indicated  by  the  cusps  in  the  contours  at  0  —  60  degrees,  which  corresponds  to  the  familiar 
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degeneracy  associated  with  D^h  symmetry.  That  is,  comparisons  of  the  surfaces  in  panels 
(a)  and  (b)  of  Fig.  4  shows  the  two  surfaces  to  pass  through  one  another,  the  former 
depicting  the  higher-value  eigensurface  to  the  right  of  the  60  degree  line  and  the  lower- 
value  eigensurface  to  the  left  of  this  line  of  intersection.  Correspondingly,  panel  (b)  shows 
the  reversed  arrangement,  with  the  lower  eigenvalue  surface  to  the  right  and  the  higher  to 
the  left  of  the  60  degree  line.  The  associated  more  familiar  energy  surfaces  for  the  ground- 
and  first-excited  H3  states  for  C-iv  atomic  arrangements  obtained  from  Eqs.  (53a)  and 
(53b)  shown  in  Fig.  5  similarly  depict  the  expected  D3h  degeneracy  and  the  associated 
surface  crossing  along  the  0  —  60  degrees  line  formed  by  the  cusps  in  the  contours,  in 
a  manner  complementary  to  the  corresponding  seam  of  intersection  anticipated  by  the 
metric-matrix  eigenvalues  of  Fig.  4. 

In  Fig.  6  are  shown  comparisons  of  the  two  doublet  energies  of  Eqs.  (53a)  and  (53b)  with 
coresponding  valence-bond  calculations  for  the  Civ  atomic  arrangements  of  Fig.  5  in  which 
Rab  =  Rac  —  1-64  a 0  are  held  fixed  and  9  is  varied  from  20  to  180  degrees.  Evidently,  the 
present  results  for  both  surfaces  are  in  good  agreement  with  the  corresponding  Is3  valence- 
bond  calculations  which  include  only  covalent  structures  (o),  and  are  also  in  accord  with 
valence-bond  results  which  include  all  possible  ionic  configurations  (•)  in  the  Is3  basis. 
The  crossing  of  the  two  surfaces  when  0  =  60  degrees  is  evident  in  the  figure,  and  in 
good  agreement  with  the  valence-bond  calculations.  As  in  Fig.  3,  the  apparent  upper 
bound  provided  by  the  present  results  on  the  valence-bond  calculations  in  the  case  of  the 
ground-state  surface  in  Fig.  6  is  a  happenstance  in  view  of  the  small  basis  representation 
employed,  with  a  similar  apparent  upper  boundedness  not  present  in  the  case  of  the  excited 
doublet  state. 

Finally,  Figs.  7  and  8  depict  contours  of  constant  S-matrix  eigenvalues  of  Eq.  (45a)  and 
energies  of  Eqs.  (53a)  and  (53b)  for  T-shaped  Civ  arrangements  in  H3  in  which  both  the 
width  (w)  and  the  height  ( h )  of  the  T  are  varied.  Evidently,  both  sets  of  these  contours 
depict  the  seams  of  intersection  along  the  line  w  —  1.115  h  of  D3/,  symmetry  associated 
with  the  surface  crossings  as  the  degenerate  configuration  is  reached,  confirming  that  the 
intersection  of  the  S-matrix  surfaces  again  anticipates  the  corresponding  crossing  in  the 
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more  familiar  energy  surfaces.  As  in  Figs.  4  and  5,  the  two  sets  of  eigenvalue  surfaces  in 
Figs.  7  and  8  depict  the  higher-  and  lower-valued  portions  of  two  smooth  surfaces  which 
intersect  at  the  indicated  seam.  Quantitative  comparisons  of  these  results  for  fixed  w  =1.64 
ao  as  a  function  of  h  with  the  valence-bond  calculations  shown  in  Fig.  9  confirms  again 
the  good  agreement  obtained  between  the  present  results  and  corresponding  valence-bond 
calculations  in  this  simple  example. 

4.0  Comparative  and  Concluding  Remarks 

The  common  features  of  previously  devised  semi-empirical  approaches  for  incorporating 
fragment  information  in  electronic  structure  calculations  are  indicated  in  Section  4.1,  de¬ 
tailed  analyses  of  the  atom-in-molecules  [13]  and  diatomics-in-molecules  [14]  methods  in 
the  context  of  the  present  development  are  reported  in  Sections  4.2  and  4.3,  respectively, 
and  concluding  remarks  in  summary  of  the  present  report  are  provided  in  Section  4.4. 

4.1  Fragment-Based  Semi-Empirical  Electronic  Structure  Approaches 

The  aims  and  objectives  of  the  present  ab  initio  development  are  closely  related  to  those 
of  earlier  semi-empirical  valence-bond-like  attempts  to  incorporate  atomic  and  diatomic 
information  in  studies  of  the  electronic  structures  of  molecules  [13,14],  The  development  of 
the  foregoing  Sections  2  and  3  provides  an  analytical  perspective  from  which  to  critically 
examine  the  theoretical  foundations  of  these  earlier  methods. 

The  theoretical  formulations  of  the  particular  semi-empirical  methods  of  interest  here  are 
based  on  adoption  of  the  ultimately  singular  Hamiltonian  matrix 

H a(R)  =  {PA&(r  :  R)\H(r  :  R)\PA<$>(r  :  R))  (54) 

constructed  in  a  non-othogonal  explicitly  antisymmetric  atomic  spectral-product  represen¬ 
tation,  and  subsequent  solution  of  the  matrix  Schrodinger  equation  [16] 

H a(R)  •  Ua(-R)  =  S (R)  •  UA (R)  ■  E (R).  (55) 

Commutation  of  the  operators  Pa  and  H(r  :  R)  in  the  Hamiltonian  matrix  of  Eq.  (54) 
gives 

Ha (R)  =  (PA<&(r  :  R)\PA\H(r  :  R)$(r  :  R)) 

This  page  is  Distribution  A:  approved  for  public  release;  distribution  unlimited.  Page  35 


=  Q1/2{PA<t>{r  :  R)\H(r  :  R)${r  :  R)),  (56) 

with  the  atoms-in-molecules  [13]  and  diatomics-in-molecules  [14]  approaches  employing 
different  approximations  to  the  vector  H[r  :  R)&(r  :  R)  in  Eq.  (56).  The  notational 
conventions  of  the  present  development  are  employed  in  Eqs.  (54)  to  (56)  in  order  to 
present  the  earlier  developments  in  forms  most  convenient  for  comparative  purposes. 

4.2  Atoms-in-Molecules  Method 

The  atoms-in-molecules  method  of  Moffitt  [13,36]  has  been  the  subject  of  early  evaluations 
[37,38]  and  continuing  but  largely  indecisive  refinements  [39-43].  In  the  original  theoretical 
formulation  [13],  the  aforementioned  vector  in  Eq.  (56)  is  given  by  the  expression  H(r  : 
R)$>(r  :  R)  =  «3>(r  :  R)  •  W  +  V(r  :  R)&(r  :  R),  where  the  matrix  W  = 
is  the  sum  of  atomic  Hamiltonian  matrices  of  Eq.  (7)  and  the  operator  V{r  :  R)  = 
)C/3=7v+i  V(a,,3\i,j  :  Rap)  is  the  sum  of  the  atomic-pair  interaction  potentials  of 
Eq.  (5).  This  expression  is  an  identity  when  the  atomic  spectral  eigenstates  employed 
are  exact  solutions  of  their  respective  atomic  Schrodinger  equations,  in  which  case  the 
associated  Hamiltonian  matrix  of  Eq.  (56)  becomes 

Ha(-R)  =  S (R)  •  W  +  Q1/2(PA$(r  :  R)\V(r  :  R)\$(r  :  R)).  (57) 

In  the  absence  of  exact  atomic  eigenstates,  a  semi-empirical  attitude  is  adopted  in  which 
Eq.  (57)  is  used  to  provide  an  approximation  to  the  Hamiltonian  matrix  of  Eq.  (54) 
by  employing  experimental  atomic  energies  in  the  diagonal  matrix  W  and  suitable  cor¬ 
responding  atomic  eigenstates,  introduced  in  an  attempt  to  provide  reliable  estimates  of 
the  Coulombic  interaction  terms  in  Eq.  (57)  which  would  be  obtained  using  the  correct 
atomic  eigenstates,  were  they  available  [13]. 

Although  the  ideas  motivating  the  atoms-in-molecules  approach  have  considerable  merit, 
the  development  of  Eqs.  (54)  to  (57)  is  flawed  in  a  number  of  fundamental  aspects  which 
have  not  been  addressed  in  the  attempted  refinements  reported  to  date  [36-43].  Specifically, 
the  two  matrices  on  the  right-hand  side  of  Eq.  (57)  are  generally  not  Hermitian,  even  in 
the  case  of  use  of  exact  atomic  spectral  eigenstates,  so  that  neither  have  any  physical 
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significance.  Moreover,  ad  hoc  enforcement  of  individual  Hermiticity  of  the  two  matrices 
results  in  a  commonly  employed  Hermit ian  expression  (l/2)(S(i?)-W+W'  -S(R) ')  for  the 
atomic  terms  which  does  not  avoid  the  presence  of  the  ultimately  singular  metric  matrix. 
Accordingly,  when  the  antisymmetric  basis  becomes  linearly  dependent,  the  solution  of 
Eq.  (55)  becomes  unstable.  Although  the  consequences  of  this  instability  have  been  often 
noted,  the  generality  and  extent  of  the  linear  dependence  of  an  antisymmetrized  atomic 
spectral-product  basis  does  not  appear  to  have  been  fully  appreciated  in  the  literature.  In 
Fig.  10  are  shown  as  an  example  the  eigenvalues  of  the  metric  matrix  of  Eq.  (10)  for  H3 
molecules  in  symmetric  collinear  arrangements  constructed  in  2slp  and  3s2p  Slater  basis 
sets  following  the  development  of  Section  3.3.  Evidently,  even  in  the  smaller  2slp  basis  set 
employed  in  panel  (a)  there  are  a  significant  number  of  eigenvalues  near  the  lower  limiting 
point  of  accumulation  of  the  spectral  interval  approriate  for  H3  (0-6),  which  clustering 
persists  even  at  larger  values  of  interatomic  separation  Rab  =  Rac.  In  the  bigger  3s2p  basis 
set  of  panel  (b),  the  S- matrix  eigenvalues  are  seen  to  become  increasingly  dense  at  the  lower 
limiting  point,  in  accordance  with  the  clustering  expected  at  the  zero  point  of  accumulation 
in  the  spectrum  of  a  compact  operator  [25].  This  common  feature  of  compact  operators 
can  be  expected  to  appear  in  the  atomic  spectral-product  representations  of  aggregate 
electron  permutation  operators  in  cases  of  many-electron  atoms  quite  generally. 

It  is  clear  that  the  metric  matrices  in  both  basis  sets  of  Fig.  10  are  highly  singular, 
requiring  elimination  of  the  linearly-dependent  basis  states  by  adoption  of  canonical  or 
other  orthogonalization  procedures.  Attempts  to  isolate  the  non-singular  portion  of  the 
Hamiltonian  matrix  of  Eq.  (57)  in  this  way,  however,  result  in  an  unfortunate  sensitivity  to 
the  atomic  energies  employed  in  the  W  matrix  appearing  there  [39-43].  Additionally,  use 
of  the  aforementioned  Hermitian  form  of  the  atomic  terms  in  Eq.  (57)  generally  renders 
largely  invalid  attempts  to  isolate  the  linearly-independent  subspace  of  the  antisymmetric 
aggregate  basis  by  unitary  transformation  [43],  suggesting  that  only  the  smallest  possible 
basis  sets  can  be  reliably  employed  in  the  atoms-in-molecules  development  of  Eqs.  (54)  to 
(57)  in  the  absence  of  modification. 

A  more  fundamental  difficulty  associated  with  the  Hamiltonian  of  Eqs.  (56)  than  that 
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indicated  in  the  foregoing  is  revealed  in  the  course  of  attempts  to  intepret  the  Moffitt 
development  of  Eqs.  (54)  to  (57)  as  a  computationally  viable  approach  in  which  the 
spectral-product  Hilbert  space  d* ( r  :  R)  potentially  has  a  limit  of  completion.  Writing  the 
vector  in  Eq.  (56)  in  this  case  in  the  form  H(r  :  R)&(r  :  R )  — »  <J>(r  :  R)  ■  H(i7),  which 
expression  is  exact  in  the  closure  limit,  gives  Ha(-R)  — * ►  S (R)  ■  H(.R),  where  H(i?)  is  the 
spectral-product  Hamiltonian  matrix  of  Eq.  (7).  In  this  case,  Eq.  (55)  becomes 

S (R)  •  H (R)  •  Ua(-R)  =  S (R)  •  Ua(-R)  ■  E(fl),  (58) 

from  which  the  metric  matrix  S (R)  might  be  factored  were  it  not  for  the  fact  that  it  is 
generally  singular  and,  consequently,  does  not  have  an  inverse. 

Since  the  matrices  S (R)  and  H(i?)  in  Eq.  (58)  generally  fail  to  commute  in  any  finite 
representation,  the  Hamiltonian  S (R)  •  H(i7)  appearing  there  is  not  Hermitian,  in  spite  of 
the  fact  that  Ha(-R)  of  Eq.  (56)  is  Hermitian  and  the  operators  Pa  and  II (r  :  R)  commute. 
That  is,  the  atomic  spectral-product  basis  does  not  provide  a  faithful  representation  of  the 
commuting  operators  Pa  and  H(r  :  R)  in  the  form  of  the  matrices  S (R)  and  H(i7), 
rendering  the  development  of  Eqs.  (54)  to  (58)  invalid  in  any  finite  representation.  This 
failure  follows  from  the  fact  that  the  operation  Pa  ( r  :  R)  is  not  closed  in  any  finite  Hilbert 
space  3>(r  :  R)  of  Eq.  (3).  In  the  closure  limit,  however,  Pa^{v  :  R)  — >  (5_1/2<J>(r  :  R)  ■ 
S (-R),  and  a  unitary  transformation  XJ(R')  can  be  found  which  simultaneously  diagonalizes 
both  the  Hamiltonian  matrix  H(i2)  of  Eq.  (7)  and  the  metric  matrix  S (R)  of  Eq.  (10), 
in  accordance  with  the  commutation  of  the  two  associated  operators.  In  this  event,  the 
matrices  S (R)  =  U (R)  ■  Sd(R)  ■  U (#)t  and  H (R)  -  U (R)  ■  E (R)  ■  U(i7)t  commute,  with 
the  spectrum  of  Sd{R)  comprising  the  lower  and  upper  end  point  vaules  0  and  Q  of  Eq. 
(9)  in  this  limit. 

In  spite  of  the  aforementioned  complications,  the  non-singular  portion  of  the  Hamiltonian 
matrix  appearing  in  Eq.  (58)  can  be  extracted  in  finite  representations,  as  well  as  in 
the  limit  of  closure,  following  the  development  of  Section  2.0.  Specifically,  isolating  the 
physical  block  of  the  Hamiltonian  matrix  of  Eq.  (56)  in  accordance  with  Eq.  (8)  and 
setting  Ha(-R)  — ►  S(R)  ■  H{R)  there  gives 

{HS(R)}PP  =  .  {Us (R)t  •  Ha(K)  ■  Us(Ji)}pp  •  {S„(fl)}„1/2 
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=  {Sd(R)}^/2  ■  {Us (R)t  ■  S(R)  ■  H (R)  ■  US(K)}„,  •  {Sd(R)}„U2 

=  {Sd(R)}+1-,2  .  {Us(fl)t  •  H(fl)  ■  Us(fl)}„  ■ 

-*{Us(H)r-H(fl)-Us(H)}„,.  (59) 

Equation  (59)  is  seen  to  constitute  a  similarity  transformation  of  the  Hamiltonian  matrix 
of  Eq.  (8)  in  finite  orders  which  converges  to  the  physical  block  of  the  spectral-product 
Hamiltonian  matrix  of  Eq.  (8)  in  the  closure  limit,  evidently  resulting  in  a  computational 
approach  essentially  identical  with  that  of  Eqs.  (1)  to  (12).  Accordingly,  it  is  seen  that 
a  strictly  Hermitian  matrix  representation  suitable  for  computational  applications  is  ob¬ 
tained  from  Mofhtt’s  development  when  its  non-singular  terms  are  isolated  employing  Eq. 
(59),  in  which  event  a  closure  limit  can  be  properly  attained  in  accordance  with  the  spec¬ 
tral  theory  of  Sections  2  and  3,  and  the  distinction  between  use  of  explicitly  antisymmetric 
and  spectral-product  basis  states  becomes  inconsequential. 

4.3  Diatomic-in-Molecules  Method 

The  diatomics-in-molecules  approach  has  been  the  object  of  continuing  applications  and 
refinements  through  the  years  [44-56].  In  the  most  commonly  employed  realizations  of  the 
original  development  of  Ellison  [14] ,  the  vector  in  Eq.  (56)  is  written  [46-50] 

H(r  :R)i>(r:R)  =  J2  E  {  3  ■  «))  }#(r  :  R), 

a=l/3=a+l  ^  ' 

(60) 

employing  Eq.  (13)  and  the  notation  of  the  present  development.  In  this  case,  Eq.  (56) 
becomes 

N  N 

H a(R)  =  Q1/2J2  E  {(PA®(r  :  R)\H^\iJ  :  Raf3)\<t>(r  :  R)) 

a=l  /3=a+ 1 

:  R)\Hfa(i)  +  H^Umr  :  H))} 

N  N  ( AT 

=  S(Jt)  ■  £  £  (61) 

a=l  /3=a+l  ' 
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where  H {^pRap)  =  (P^,/3)$(r-  :  R)\H^\i,  j  :  R^P^  $(r  :  R))  =  (Q(Q^)1/2 
(Pf’^#(r  :  R)\H^(i,j  :  Rap)\<t>(r  :  R))  is  the  diatomic  Hamiltonian  matrix  in  the 
lion-orthogonal  antisymmetrized  spectral  product  pair  basis  P'|k’/5'<h(r  :  R) .  The  atomic 
integrals  in  the  first  line  of  Eq.  (61)  have  been  evaluated  by  employing  the  expression 
(H^pi)  +  ffW(j))$(r  ;  R)  =  $(r  :  R)  •  (H^  +  H^),  whereas  the  diatomic  integrals 
are  evaluated  in  the  form 


Q1/2(PA®{r  :  R)\H^\i,j  :  Rap)\<t>(r  :  R)) 

=  Q1/2(Q^)-1/2(PA^(r  :  R)\H^\i,j  :  Rap)\P^P)^{r  :  R)) 

=  Q1/2(Q{a’p))~1/2{PA$(r  :  R)\P{°’®*(r  :  B)>  •  S^J^)"1  •  H {^P)(Rap) 

=  S (R)  ■  S^iR^)-1  ■  n^iRap)  (62) 

using  the  representation  of  H^a,^pi,j  :  Rap)P^‘,^^(r  :  R)  in  the  basis  P^*’^<I>(r  :  i?)  in 
the  form  H^piJ  :  JRQ/3)pJ’/3)$(r  :  i?)  -  pj,/3)#(r  :  i?)  •  S^^a/j)-1  -H^pRap). 

The  Schrodinger  equation  for  the  diatomic  Hamiltonian  matrix  H^’^(Hayg), 

H^i*^)  ■  -  S^\Rap)  •  U {^pRap)  •  E^^(i2Q/3),  (63) 

can  be  used  to  formally  eliminate  the  matrix  S^^pRap)-1  in  Eq.  (62)  by  writing 
S^pRap)-1  ■  U^pRap)  =  U i'P){Rap)  •  E ^pRap)  ■  U^pRap)-1,  in  which  form 
Eq.  (62)  is  commonly  employed  [46-50].  The  matrices  S ^a’^pRap)  and  E ^a^pRa  d)  of 
Eqs.  (61)  to  (63)  depend  explicitly  on  both  the  magnitude  and  direction  of  the  separation 
vector  Rap,  as  do  the  matrices  H^’^pRap)  and  (Rap)  of  Eq.  (63),  in  contrast 

to  their  counterparts  of  Eqs.  (24)  and  (26),  which  are  defined  in  terms  of  the  rotated 
representations  of  Eq.  (19). 

The  development  of  Eqs.  (60)  to  (63)  entails  commutation  of  both  aggregate  and  di¬ 
atomic  Hamiltonian  operators  with  the  corresponding  antisymmetrizers,  resulting  in  the 
appearance  of  both  aggregate  and  diatomic  metric  matrices  and  the  non-Hermitian  forms 
of  the  atomic  and  diatomics  Hamiltonian  matrices  in  Eq.  (61),  the  latter  consequent  of  the 
failure  of  the  metric  and  Hamiltonian  matrices  to  commute  in  any  finite  representation. 
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Accordingly,  instabilities  associated  with  the  linear  dependencies  arising  from  the  aggre¬ 
gate  and  diatomic  metric  matrices  can  be  expected  to  appear  in  applications  based  on  the 
non-Hermitian  diatomics-in-molecules  expression  of  Eq.  (61).  Specifically,  a  sensitivity  to 
input  atomic  and  diatomic  information  indicated  in  connection  with  Eqs.  (57)  and  (58) 
above  also  applies  to  the  diatomics-in-molecules  Hamiltonian  matrix  of  Eqs.  (60)  to  (63). 
Attempts  to  deal  with  the  consequences  of  these  circumstances  include  the  widespread 
avoidance  of  the  aggregate  metric  matrix  by  the  ad  hoc  approximation  S (R)  — >  I,  adoption 
of  explicitly  Hermitized  forms  of  the  atomic  and  interaction  matrices,  and  dealing  directly 
with  the  consequences  of  a  non-Hermitian  Hamiltonian  formalism.  Although  many  of  the 
applications  of  Eq.  (61)  have  shown  considerable  inventiveness  in  overcoming  computa¬ 
tional  obstacles  [44-56],  the  inherently  flawed  nature  of  the  development  has  apparently 
precluded  widespread  adoption  of  the  approach. 


As  in  the  case  of  the  development  of  Section  4.2,  a  computationally  viable  diatomics- 
in-molecules  approach  can  be  based  on  the  linearly-independent  subspaces  of  explicity 
antisymmetric  aggregate  and  diatomic  spectral-product  representations  of  Eqs.  (40)  and 
(41).  Specifically,  the  non-singular  block  of  Eq.  (61)  is  isolated  in  the  form  [cf. ,  Eq.  (59)] 

{HS(JJ)}PP  =  {Sd(fi)}-y2  .  {Us(R)t  .  HA(Jl)  ■  Us(fl)}„  ■  {Sd(H)}W2 

N  N 

=  {Sd(«)}+1/2  ■  XI  E 

cx.=l  /3=a+l 


Ka)(fl)  +  .  {Sd(ii)},-1/2 


(64) 


I  PP  ^  ~  '  'JPP  ’ 

~  / 

where  HgV^(i7)  and  Hg f\R)  are  the  transformed  atomic-energy  matrices  of  Eq.  (8)  and 
the  subscript  “pp”  on  the  diatomic  terms  refers  to  both  the  aggregate  and  diatomic  rep¬ 
resentations.  The  latter  terms  in  Eq.  (64)  can  be  evaluated  by  inserting  the  identities 
uf’/3)(JRa/3)-Uf’/3)(JRa^)t  =  V{sa^(Ra^-V{sa’p\RaP)  =  I  before  after  S ^\Rap)~\ 
respectively,  and  noting  that 

S ^{Rap)-1  =  uf’/3)(i?ap)t  •  S •  V^\Raf3)  (65) 

when  restricted  to  the  indicated  linearly-independent  diatomic  subspace,  to  obtain 

S ^{Rafi)-1  ■  H £,I3\R*I3)  =  U {s^(Rap)  ■  S ^(Rap)-1  •  U ^  (Ra^  ■  H i^iRap) 
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=  U ^(R^)  ■  S<“’'3)(Ra|3)-1  ■  ■  H (Z’m(Ra0)  ■  U ^(Raf,)  ■  U ^(R^)' 

=  V(sA(Rat 3)  •  {S^'ffl^)}^2  ■  {H<,“"«(Ji^)}pp  •  {S {^\Ro,e)}l'p  ■ 

(66) 

Here, 

is  the  linearly-independent  block  of  the  Hamiltonian  matrix  H ^,'8',(i2Q/g),  and  only  the 
non-singular  portion  of  the  diagonal  matrix  has  been  retained  in  the  last 

line  of  Eq.  (66)  in  accordance  with  the  specification  of  Eq.  (64).  Finally,  employing  Eq. 
(66)  in  Eq.  (64)  and  taking  the  closure  limit  as  in  Eqs.  (40)  to  (42)  gives 


N  N 


{lls(R)}PP  =  E  E  +  Hf  (fl)}w}  (68) 


(-^  2)  rxj(a) 

\±i')  Spp - 

a=l  /3=a+l 

for  the  diatomics-in-molecules  Hamiltonian  matrix,  where 


{H<“-»(fl)}pp  =  {U^“'w(fl)}T  ■  {H^'(R„3)}pp  •  {Ug“’p,(J?)} 


—  rTi(a’d)/ 


it 


(a,/3) , 


r(«,/3)/ 


PP 


(69) 


and 

{l4“’»(JJ)}pp  =  {U^fi^)'  ■  Us(fl)}pp.  (70) 

Equations  (68),  (69),  and  (70)  are  seen  to  be  identical  with  Eqs.  (32),  (33),  and  (42), 
respectively,  when  note  is  taken  of  the  fact  that  the  spatial  frame  transformation  of  Eqs. 
(19)  to  (21)  has  not  been  incorporated  in  the  devlopment-  of  Eqs.  (64)  to  (70).  When  this 
transformation  is  introduced,  the  modified  diatomics-in-molecules  Hamiltonian  of  Eqs. 
(68)  to  (70)  is  seen  to  be  identical  with  that  obtained  from  the  computational  implemen¬ 
tation  of  the  spectral-theory  in  Section  3  employing  explicitly  antisymmetric  basis  states. 
Accordingly,  a  strictly  Hermitian  diatomic-in-molecules  formalism  suitable  for  computa¬ 
tional  applications  is  obtained  from  the  original  largely  semi-empirical  development  when 
the  insights  provided  by  the  spectral-theory  of  Sections  2  and  3  are  incorporated  into  the 
analysis. 
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4.4  Concluding  Remarks 


The  theoretical  development  reported  here,  in  which  antisymmetry  requirements  are  ap¬ 
plied  subsequent  to  construction  of  the  many-electron  Hamiltonian  matrix  in  an  orthogonal 
spectral-product  basis,  provides  an  alternative  perspective  on  electronic  structure  calcula¬ 
tions  of  atoms,  molecules,  and  other  forms  of  matter.  Permutation-symmetry  basis  states 
are  seen  to  play  a  central  role  in  the  development,  emphasizing  the  importance  of  the  eigen- 
spectrum  of  the  matrix  rcpresentatative  of  the  aggregate  electron  antisymmetrizer,  which 
is  seen  to  converge  to  the  two-point  spectrum  commonly  associated  with  a  compact  opera¬ 
tor.  The  approach  provides  for  the  definition  of  Hermitian  matrix  representatives  of  atomic 
and  diatomic  operators  which  individually  have  appropriate  asymptotic  separation  limits, 
and  which  accommodate  the  incorporation  of  diatomic  fragment  information  in  forms  free 
of  ad  hoc  assumptions  or  approximations.  The  spectrum  of  the  aggregate  Hamiltonian 
matrix  obtained  in  this  way  satisfies  the  Hylleraas-Undheim  theorem  for  fixed  dimensions 
of  the  permutation- symmetry-blocked  diatomic  Hamiltonian  matrices  as  the  dimensions  of 
the  aggregate  permutation-symmetry  basis  is  enlarged  in  a  linear  variational  development. 
Accordingly,  convergence  to  specific  portions  of  the  correct  Hamiltonian  eigenspectrum  can 
be  anticipated  from  the  development  employing  converged  representations  of  the  underly¬ 
ing  bare  diatomic  states  constructed  in  sufficiently  large  spectral-product  representations. 

A  particular  computational  implementation  of  the  general  development  makes  use  of  largely 
conventional  methods  of  calculations  for  atoms  and  diatomic  molecules,  and  can  take  ad¬ 
vantage  of  the  currently  available  highly  developed  methods  devised  for  this  purpose.  The 
simple  H3  example  studied  illustrates  the  important  role  of  the  eigensurfaces  of  the  metric 
matrix  in  identifying  non-Pauli  states  in  the  aggregate  spectral-product  basis,  on  the  one 
hand,  and  of  the  presence  of  linear  dependence  in  the  corresponding  antisymmetrized  basis, 
on  the  other.  Additionally,  the  S-matrix  eigensurfaces  anticipate  the  general  topological 
features  of  the  corresponding  more  familiar  potential  energy  surfaces,  and,  in  particular,  in 
case  of  the  low-lying  doublet  states  of  H3  are  seen  to  anticipate  the  presence  of  well-known 
electronic-energy  surface  crossings  at  high-symmetry  configurations.  Additional  calcula¬ 
tions  on  other  molecules  which  include  many-electron  atoms  are  required  to  confirm  the 
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promise  of  these  illustrative  applications  to  the  H3  molecule. 

The  development  further  provides  a  theoretical  basis  from  which  to  analyze  previously 
reported  highly  promising  but  largely  unrealized  methods  for  incorporating  fragment  in¬ 
formation  in  electronic  structure  determinations  [13,14].  The  failure  to  obtain  strictly 
faithful  matrix  representatives  of  the  aggregate  electron  antisymmetrizer  and  Hamiltonian 
in  the  spectral-product  basis,  consequent  of  the  absence  of  closure  in  the  operation  of  the 
antisymmetrizer  in  any  finite  spectra-product  representation,  is  seen  to  be  the  root  cause 
of  limitations  and  aberrant  behaviors  observed  in  applications  of  these  approaches  [44-50]. 
Specifically,  the  failure  of  the  matrix  representatives  of  the  antisymmetrizer  and  Hamilto¬ 
nian  operators  in  the  specral-product  basis  to  faithfully  represent  the  commutation  of  the 
defining  operators  leads  to  the  appearance  of  non-Hermitian  matrix  representatives  and 
instabilities  associated  with  the  singular  nature  of  both  aggregate  and  diatomic  metric 
matrices.  These  limiting  attributes  of  the  earlier  developments  are  overcome  employ¬ 
ing  the  more  complete  understanding  of  the  role  of  the  metric  matrices  provided  by  the 
spectral-theory  analysis.  Extraction  of  the  non-singular  portions  of  aggregate  Hamilto¬ 
nian  matrices  constructed  in  linearly-depcndent,  explicitly  antisymmetrized  basis  sets  are 
seen  to  result  in  commutationally  viable  expressions  from  both  atoms-in-molecules  and 
diatomics-in-molecules  approaches  which  are  identical  with  the  results  obtained  here,  pro¬ 
viding  confidence  in  applications  of  the  modified  versions  of  these  earlier  methods  and  of 
the  present  development. 
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Figure  1-Constant-value  contours  for  the  two  doublet  eigenvalues  (a)  {Sd(R)}i  and  (b) 
{Sd(R)}2  of  the  S(-R)  matrix  of  Eqs.  (44)  and  (45)  for  collinear  atomic  configurations 
of  H3.  The  coordinate  axes  give  the  distances  Rat,  and  Rac  between  the  central  atom 
located  at  the  origin  and  the  outer  atoms,  measured  in  atomic  units  (ao),  employing  an 
increment/decrement  of  0.1  between  adjacent  contours.  In  panel  (a)  the  contours  increase 
uniformly  in  value  from  that  labeled  1.05,  except  as  R,,},  =  Rac  — >  0  inside  the  saddle 
point,  in  which  region  the  values  decrease  monotonically  to  zero,  whereas  in  panel  (b)  the 
contour  values  decease  uniformly  from  that  labeled  0.95,  as  is  discussed  further  in  the  text. 
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Figure  2-Const  ant-energy  contours  for  the  ground  (a)  and  first-excited  (b)  doublet- 
state  potential  energy  surfaces  of  H3  for  noil-symmetric  collinear  atomic  arrangements, 
obtained  from  the  matrix  elements  {Hg(i2)}n  and  {Hs(-R)}22  of  Eqs.  (51)  and  (52). 
The  coordinate  axes  are  as  in  Fig.  1,  measured  in  atomic  units  (ao),  the  zero  of  energy 
employed  (-1.50  au)  is  that  of  the  three-atom  dissociation  limit  (H+H+H),  and  an  energy 
decrement /increment  of  0.5  eV  is  employed  between  contours.  Panel  (a)  depicts  the  fa¬ 
miliar  saddle  and  diatomic  structures  of  the  ground-state  energy  surface  in  H3,  whereas 
the  excited-state  energy  surface  of  panel  (b)  is  montonically  increasing  with  decreasing 
separations  Ra},  and/or  Rac.  as  is  discussed  further  in  the  text. 
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Figure  3-Ground  and  first-excited  doublet  state  energies  in  H3  for  symmetric  collinear 
atomic  arrangements.  The  solid  curves  depict  the  values  obtained  from  Eqs.  (53a)  and 
(53b),  respectively,  of  the  present  development,  whereas  the  data  points  are  obtained  from 
conventional  valence-bond  calculations  in  a  Is3  basis  set  including  (•)  and  excluding  (o) 
ionic  configurations.  The  zero  of  energy  (-1.50  au)  is  that  of  the  three-atom  dissociation 
limit  (H+H+H),  as  in  Fig.  2. 
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Figure  4-Constant-value  contours  for  the  two  doublet  eigenvalues  (a)  {Sd(i?)}i  and 
(b)  {S d{R)}2  of  the  S (R)  matrix  of  Eqs.  (44)  and  (45)  for  C^v  atomic  arrangements 
(Rab  —  Rac )j  where  6  is  the  apex  angle  of  the  triatomic  H3  configuration,  employing  an 
increment/decrement  of  0.1  between  adjacent  contours.  The  cusps  in  the  contours  of  both 
panels  indicate  the  presence  of  a  seam  of  intersection  in  the  two  surfaces  at  the  0  =  60 
degree  line,  with  the  contour  values  decreasing  uniformly  from  right  to  left  in  panel  (a) 
and  decreasing  uniformly  from  top  to  bottom  in  panel  (b),  as  is  discussed  further  in  the 
text. 
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Figure  5-Constant-energy  contours  for  the  two  lowest-lying  doublet  states  of  H3  for  C2,, 
atomic  arrangements  (Rab  —  Rac ),  obtained  from  Eqs.  (53a)  and  (53b),  where  6  is  the 
apex  angle  of  the  tri atomic  configuration,  employing  an  energy  increment/decrement  of 
0.5  eV  between  adjacent  contours.  The  zero  of  energy  (-1.50  au )  is  that  of  the  three-atom 
dissociation  limit  (H+H+H).  The  cusps  in  the  contours  of  both  panels  indicate  the  familiar 
seam  of  intersection  in  the  two  energy  surfaces  at  the  9  =  60  degree  line  in  Dy,  symmetry, 
with  the  contour  energy  values  increasing  uniformly  from  top  to  bottom  in  both  panels. 
Panel  (a)  depicts  the  lowest-lying  surface  only  for  large  values  of  6 ,  whereas  panel  (b) 
depicts  the  lowest-lying  surface  only  for  small  0,  as  is  discussed  further  in  the  text. 
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Figure  6-Ground  and  first-excited  doublet-state  energies  in  H3  for  C^v  atomic  arrange¬ 
ments  (Rab  —  Rac  —  1-64  ao),  where  6  is  the  apex  angle  of  the  triatomic  configuration. 
The  solid  curves  depict  the  values  obtained  from  Eqs.  (53a)  and  (53b)  of  the  present  devel¬ 
opment,  whereas  the  data  points  are  obtained  from  conventional  valence-bond  calculations 
in  a  Is3  basis  set  including  (•)  and  excluding  (o)  ionic  configurations.  The  zero  of  energy 
(-1.50  au )  is  that  of  the  three-atom  dissociation  limit  (H+H+H),  as  in  Fig.  5.  The  crossing 
point  of  the  two  curves  corresponds  to  D3h  symmetry,  as  discussed  in  the  text. 
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Figure  7-Constant-value  contours  for  the  two  doublet  eigenvalues  (a)  {Sd(R)}i  and  (b) 
{Sd(R)}2  of  the  S (R)  matrix  of  Eq.  (44)  and  (45)  for  T-shaped  atomic  arrangements 
(width  -  w,  height  -  h)  in  H3,  employing  an  increment/decrement  of  0.1  between  adjacent 
contours.  In  panel  (a)  the  contours  are  monotonically  increasing  from  that  labelled  1.0, 
whereas  in  panel  (b)  they  are  monotonically  increasing  from  that  labelled  0.9.  The  cusps 
in  the  contours  of  both  panels  indicate  the  presence  of  a  seam  of  intersection  in  the  two 
surfaces  along  the  w  —  1.155  h  line  corresponding  of  D^h  symmetry,  as  is  discussed  further 
in  the  text. 
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Figure  8-Constant-energy  contours  for  the  two  lowest-lying  doublet  states  of  H3  obtained 
from  Eqs.  (53a)  and  (53b)  for  T-shaped  atomic  arrangements  (width  -  w,  height  -  h). 
employing  an  energy  increment/decrement  of  0.5  eV  between  adjacent  contours,  with  the 
zero  of  energy  (-1.50  au )  employed  that  of  the  three-atom  dissociation  limit  (H+H+H). 
The  cusps  in  the  contours  of  both  panels  indicate  the  presence  of  a  seam  of  intersection 
in  the  two  surfaces  along  the  w  —  1.155  h  line  corresponding  to  /T,/,  symmetry.  Panel  (a) 
depicts  the  lowest-lying  surface  only  for  large  values  of  height  h,  whereas  panel  (b)  depicts 
the  lowest-lying  surface  only  for  small  h,  as  is  discussed  further  in  the  text. 
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Figure  9-Ground  and  first-excited  doublet-state  energies  in  H3  for  T-shaped  C^v  atomic 
arrangements  (width  -  w  =  1.64  a 0,  height  -  h ).  The  solid  curves  depict  the  values  obtained 
from  Eqs.  (53a)  and  (53b)  of  the  present  development,  whereas  the  data  points  are  obtained 
from  conventional  valence-bond  calculations  in  a  Is3  basis  set  including  (•)  and  excluding 
(o)  ionic  configurations.  The  zero  of  energy  (-1.50  au)  is  that  of  the  three-atom  dissociation 
limit  (H+H+H),  as  in  Fig.  8. 
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Figure  10-Eigenvalue  spectra  of  the  S-matrix  for  symmetric  collinear  H3  constructed 
in  (a)  2slp  and  (b)  3s2p  basis  sets.  The  limiting  points  of  the  spectra  are  seen  to  be 
zero  and  Q  =  ^/(ni^In^!)  —  6,  with  the  latter  eigenvalues  corresponding  to  the  totally 
antisymmetric  states  and  the  former  to  the  unphysical  (non-Pauli)  states  of  Eqs.  (9)  to 
(11).  The  lower  end  of  the  spectrum  corresponds  to  a  point  of  accumulation,  as  discussed 
in  the  text. 


